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ABSTRACT 

The tidal stirring model posits the formation of dwarf spheroidal galaxies (dSphs) via the tidal interactions 
between late-type, rotationally-supported dwarfs and Milky Way-sized host galaxies. Using a comprehensive 
set of collisionless A^-body simulations, we investigate the efficiency of the tidal stirring mechanism for the 
origin of dSphs. In particular, we examine the degree to which the tidal field of the primary galaxy affects the 
sizes, masses, shapes, and kinematics of the disky dwarfs for a range of dwarf orbital and structural param- 
eters. Our study is the first to employ self-consistent, equilibrium models for the progenitor dwarf galaxies 
constructed from a composite distribution function and consisting of exponential stellar disks embedded in 
massive, cosmologically-motivated dark matter halos. Exploring a wide variety of dwarf orbital configurations 
and initial structures, we demonstrate that in the majority of cases the disky dwarfs experience significant mass 
loss and their stellar distributions undergo a dramatic morphological, as well as dynamical, transformation. 
Specifically, the stellar components evolve from disks to bars and finally to pressure-supported, spheroidal sys- 
tems with kinematic and structural properties akin to those of the classic dSphs in the Local Group (LG) and 
similar environments. The self-consistency of the adopted dwarf models is crucial for confirming this com- 
plex transformation process via tidally-induced dynamical instabilities and impulsive tidal heating of the stellar 
distribution. Our results suggest that such tidal transformations should be common occurrences within the cur- 
rently favored cosmological paradigm and highlight the key factor responsible for an effective metamorphosis 
to be the strength of the tidal shocks at the pericenters of the orbit. We also demonstrate that the combination of 
short orbital times and small pericentric distances, characteristic of dwarfs being accreted by their hosts at high 
redshift, induces the strongest and most complete transformations. Our models also indicate that the efficiency 
of the transformation via tidal stirring is affected significantly by the structure of the progenitor disky dwarfs. 
While the mass-to-light ratios, M/L, of the dwarf galaxies typically decrease monotonically with time as the 
extended dark matter halos are efficiently tidally stripped, we identify a few cases where this trend is reversed 
later in the evolution when stellar mass loss becomes more effective. We also find that the dwarf remnants 
satisfy the relation V max = v^cr*, where cr* is the one-dimensional, central stellar velocity dispersion and V max 
is the maximum halo circular velocity, which has intriguing implications for the missing satellites problem. 
Assuming that the distant dSphs in the LG, such as Leo I, Tucana, and Cetus are the products of tidal stirring, 
our findings suggest that these galaxies should have only been partially stirred by the tidal field of their hosts. 
We thus predict that these remote dwarfs should exhibit higher values of V mt /a*, where V lot is the stellar rota- 
tional velocity, compared to those of dSphs located closer to the primary galaxies. Overall, we conclude that 
the action of tidal forces from the hosts constitutes a crucial evolutionary mechanism for shaping the nature 
of dwarf galaxies in environments such as that of the LG. Environmental processes of this type should thus be 
included as ingredients in models of dwarf galaxy formation and evolution. 

Subject headings: galaxies: dwarfs - galaxies: fundamental parameters - galaxies: kinematics and dynamics - 
galaxies: Local Group - galaxies: structure - cosmology: dark matter 



1. INTRODUCTION 

The currently favored cold dark matt er (CDM) paradigm 
of hierarchical structur e formation (e.g.. lWhite & Reeslll978l : 
iBlumenthal et alj[l984l) generically predicts that dwarf galax- 
ies comprise the primary building blocks of more mas- 
sive systems. Observational evidence has recently con- 
firmed this expectation with the discovery of streams 
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and complex stellar structures associated with accreted 
and tidally disrupted dwarfs in the Milky Way (MW) 
(e.g.. llbata etalJ[l994l lYannv et alj|2000t llbata et aLlfe OOlb; 
Newb erget all I2002T iMaiewski et all 12003b iMarfin et all 
2004t iMartmez-Delgado et alj 120051: iGrillmair & Dionatosl 
1200a iBelokurov et all 120061). the Andromeda galaxy (M31) 
(llbata et alj 12001 at iFerguson eTaT] 120021 l2005HKalirai et all 
2006; Ibat a et all 120071), and beyond the Local Group (LG) 
(e.g. [Malin & Hadlevl Il997t IShang et all 119981: iPeng et alj 
120021: Forbes et alj|2003h iPohlen et al.ll2004 r~Understanding 
the formation and evolution of dwarf galaxies is crucial for 
testing the predictions of the CDM theory and gaining insight 
into the physical processes of structure formation. 

In this context, the dwa rf spheroidal galaxie s (dSphs) of the 
LG (see lMatedfl998l and iTolstov et al.ll2009l for comprehen- 
sive reviews) constitute excellent candidates for constraining 
the CDM model, as they are believed to be highly dark matter 
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(DM) dominated, with mass-to-light ratios of M/Ly ~ 10 1 3 
(e.g..lMateolll99aiGilmore et alJl2007HSimon & Gehall2007l 



IWalker et al. 120091) (for alter native explanations o f their ex- 
treme DM content see, e.g., lKuhn& Milled [T989I: Hlgrom 
[1991 iKroupal [19971 iLokai 12001 7" dSphs are the faintest 
galaxies known and their stellar components are supported by 
random motions, with a ratio of rotational velocity t o line- 
of-sig ht, central velocity dispersion of V m t/a t < 1 dMateol 
1998). Among the dwarf galaxies in the LG, dSphs are also 
the most numerous and they tend to be clustered around the 
massive spirals MW and M3 1 (though some outliers exist in- 
cluding the distant dwarf Leo I and the isolated dSphs Ce- 
tus and Tucana that lie on the outskirts of the LG). This 
tendency is referred to as the morphology-density relation 
where rotationally-supported dwarf irregular g alaxies (dlrrs) 
are found in the periphery of the LG (e.g., Mateo 1998; 
lGrebellfl999h . Mor eover, dSphs are gas poor or completely 
devoid of gas (e.g., lMateolll998l; iGrcevich & Putmanll2009l) 
and they ex hibit a wide diversity in their st ar formation his- 
tories (e.g.. lGrebelll2000t lOrban et al.ll2008l) . Owing to their 
proximity, dSphs have been most thoroughly studied in the 
LG. However, dwarfs with the properties of dSphs have been 
rec ently identified in nearby groups and cl usters of galax- 
ies dChiboucas et alj|2009t iPennv et al1 l2009) suggesting that 
these intriguing galaxies are not unique to the LG. 

While our understanding of dSphs has grown impres- 
sively in the past decade, a de finitive mode l for their ori- 
gin s till remains elusive (see iMaverl 1201 Ol and iRravtsovl 
120101 for recent reviews). Two main classes of models 
have been proposed so far to explain their formation and 
present-day properties. In the first, dSphs are the result of 
the interplay between cosmic reonization and stellar feed- 
back suppressing g as accretion and star formation in low- 
mass galaxies (e.g.. iDekel & Silklll986t iBullock et alj|2000t 
Susa & Umemurall2004t IRicotti & Gnedinll2005t ITassis et all 
20081: ISawala et al.1 120101: iMaccio et all 120101) . In the sec- 
ond, the origin of dSphs is intimately linked to vari- 
ous environmental mech anisms, including t i dal and ram 
pressure stripping (e.g., Einasto et al.l 1 19741: iFaber & Lin ' 



19831: iMavereta l. 2001ab; Kr avtsov e^j2004lMavCTetal] 
20061 120071 IKlimentowski et al I 120071: iPefiarrubia et al.l 
2008; Klimentowsk i et all l2009al) and resonant stripping 
dD'Onghia et alj|2009t) . 

Using controlled simulations of individual disky dwarf 
galaxie s orbiting inside a MW-sized host, iMaver et al.l 
(200 lab) demonstrated for the first time that the repeated 
action of tidal forces from the primary galaxy can trans- 
form the dwarfs into pressure-supported stellar systems with 
the structural and kinematic pro perties of dSphs (see also 
IKlimentowski et aTH 2007, 2009a). This transformation mech- 
anism, termed "tidal stirring", involves a combination of 
tidally-induced dynamical instabilities in stellar disks (e.g., 
the bar and buckling instabilities) and impulsive tidal heat- 
ing of the stellar distribution. In the context of CDM, 
tidal stirring should be particularly effective since satel- 
lites of massive galaxies are affected by strong tidal forces 
due to their highly e ccentric orbits (e.g., iGhigna et al.lll998t. 
iDiemand et al.ll2007l) . 



Although the tidal stirring model naturally explains the ten- 
dency of dSphs to be concentrated near the dominant spiral 
galaxies, it is only applicable to classic dSphs, namely those 
that were known before the discovery of the ul tra-faint dwarfs 
by the Sloan Digital Sky Survey (SDSS) (e.g JSimon & Gehal 
20071). When combined with ram pressure stripping and the 



effect of radiation fields at high redshift, such as the cosmic 
ionizing ultraviolet background, the tidal stirring mechanism 
can also account for both the low gas fraction and the ex- 
tremely high DM conten t in some of the classic d Sphs such 
as Draco and Ursa Minor ( IMaver et al.l2006l [2007b . Recently, 
the structure and kinematics of the stellar core of the nearby 
Sagittarius dwarf galaxy have also bee n successfully mode led 
within the framework of tidal stirring (Lokas et al. 20 10bl) . 

In addition to controlled numerical experiments, signif- 
icant theoretical effort has been devoted to performing 
fully cosmological, hyd rodynamical simulations of dwarf 
galaxy formation (e.g., IRicotti & G nedin 2005; Rea d et al 
2006a; ITassis et al.l20 08; Governato et al. 2010; Sawa la et al 



I2010B . While pressure-supported objects with the proper- 
ties of LG dSphs have been produced in so me cases, (e.g., 
IRicotti & Gnedinl 12005b ISawala et all 120101) . other studies 
have demonstrated the formation of systems with properties 
similar to those of the likely progenitors of dSphs a ccord- 
ing to the tidal stirring model (iGovernato et al.1 12010). Suc- 
cesses notwithstanding and despite the continuing increase in 
dynamic range, limited resolution prohibits current cosmolog- 
ical simulations with hydrodynamics to address in detail the 
dynamical and structural evolution of dwarf satellite galaxies. 
In addition, due to numerical loss of angular momentum and 
overcooling, satellites in cosmological simul ations of galaxy 
forma tion are too bright and too dense (e.g., IGovernato et al.1 
l2007h . 

Given these facts, we are motivated to explore the tidal 
evolution of individual, rotationally-supported dwarf galaxies 
inside MW-sized hosts using a large ensemble of controlled 
Af-body simulations. While our work is informed by many 
past numerical investigations aimed to elucidate the tidal st ir- 
ring scenario for the origin of dSphs (IMaver et al.1 12001 alibi 
2002; Kli mentowskfetllll 120071 12009a). our numerical ex- 
periments extend those of earlier studies in several impor- 
tant respects. For example, unlike previous work, we employ 
self-consistent, equilibrium numerical models of disky dwarfs 
for our experiments. These models are derived from three- 
integral composite distribution functions (DFs) and are thus 
superior to those constructed from simpler approximate pre- 
scriptions tha t produce model s that are initially out of equi- 
librium (e.g., iHernquisllI 1993b . The transformation via tidal 
stirring is an intricate process with several distinct stages that 
depends sensitively on the development of tidally-induced dy- 
namical instabilities (e.g., bars) and the detailed response of 
the stellar distribution to tidal shocks. For this reason, the 
self-consistency of the adopted models is crucial. 

Moreover, we explore a wide variety of orbital configura- 
tions and structural parameters for the progenitor dwarf galax- 
ies, conducting a simulation campaign that allows the inves- 
tigation of a much larger parameter space than before. Our 
ultimate goal is to determine the degree to which the evolu- 
tion of rotationally-supported dwarfs is affected by the strong 
tidal field of their host galaxies under a broad range of ini- 
tial conditions, and through this to establish the generic effi- 
ciency of the tidal stirring mechanism for the origin of dSphs. 
Lastly, our numerical experiments are characterized by much 
higher numerical resolution than that of previous related stud- 
ies. This fact, in conjunction wit h the increasing a ccuracy 
of current observational data (e.g., Walker et al. 2009), offers 
unique opportunities for a systematic and quantitative com- 
parison with observations, and we undertake such a task in a 
companion paper (Lokas et al. 2010, in preparation). 

In this study, we focus on the evolution of the intrinsic, 
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global parameters of accreted disky dwarfs orbiting inside a 
MW-sized primary galaxy. Our results firmly establish that 
tidal interactions between late-type dwarf galaxies and their 
hosts can produce objects with kinematic and structural prop- 
erties akin to those of the classic dSphs in the LG and sim- 
ilar environments under a wide variety of initial conditions. 
We conclude that such tidal transformations should be com- 
mon occurrences within the currently favored cosmological 
paradigm. Environmental mechanisms of the type highlighted 
in the present work should thus constitute important ingredi- 
ents in models of dwarf galaxy formation and evolution. 

The outline of this paper is as follows. In §|2] we introduce 
the dwarf and host galaxy models. In this section, we also 
describe the methods adopted and the setup of the numerical 
experiments performed in the present study. In §[3]and §|4] we 
investigate the degree to which the orbital parameters and ini- 
tial structure of the progenitor disky dwarfs can influence their 
morphological and dynamical transformation into dSphs. Im- 
plications and extensions of our findings are presented in §|5J 
which also includes a discussion of the caveats of the current 
study and a number of promising directions for future work. 
Lastly, in §|6]we summarize our main results and conclusions. 

2. SIMULATIONS AND METHODS 

2.1. Models of Dwarf Galaxies 

The present investigation utilizes fully self-consistent, equi- 
librium models of dwarf galaxies for tidal stirring experi- 
ments. This aspect of the modeling constitutes one of the 
major improvements we introduce in this study. As we dis- 
cussed in the previous section, the nature of the transfor- 
mation mechanism is such that the self-consistency of the 
adopted dwarf mo dels is essential. Specifica lly, we em- 
ploy the method of Wid row & Dubinskf d2005l) to construct 
numerical realizations of dwarf galaxies consisting of ex- 
ponential s tellar disks embedd ed in cuspy, cosmologically- 
motiv ated iNavarro et alj d 19961 hereafter NFW) DM halos. 
The IWidrow & Dubinski d2005l) models are specified by a 
large number of parameters. They are derived from three- 
integral, composite DFs and thus represent self-consistent, 
equilibrium solutions to the co upled Poisson and co l lision- 
less Boltzmann equations. The IWidrow & DubinskH (l2005h 
method has been recently used in a variety of numerical 
studies associated with instabilities in disk galaxies, includ- 
ing the dynamics of warps and bars ( Du binski & Chakrabartvl 
l2009tlDubinski et al.ll2009h and the heating of galactic disks 
by halo substructure dGauthier et all 120 06: Kazantzid is et alj 
120081: iPurcell et alj|2009t iKazantzidis et alj|2009f) . We refer 
the reader to IWidrow & Dubin ski ( 2005) for an overview of 
all relevant parameters and a detailed description of this tech- 
nique. 

The density distributions of the dwarf DM halos are given 
by 

f>NFw(r)=- — — -J, (1) 

( r /r s )(l + r /r s ) 2 

where p s is a characteristic inner density, and r s denotes the 
scale radius of the density profile defined as the distance from 
the center where the logarithmic slope, din p(r) /dlnr, is equal 
to -2. The NFW density profile is formally infinite in extent 
with a cumulative mass t hat diverges as r — > oo. In order to 
keep the total mass finite, Wid row & DubinskH (12005) impose 
a tidal radius in the DM halo, /?/,, which represents the outer 
edge of the system. In our modeling, R/, is chosen in such 



a way that it becomes roughly equal to the cosmologically- 
motivated virial radius, R/ t ps r v ; r . As a result of this choice, 
the total mass of the halo within /?/„ which we denote M/,, 
would be equivalent to the virial mass. We control the shape 
of the halo density profile via the concentration parameter c = 
fva/r s ps R/,/r s . Higher values of concentration correspond to 
a larger fraction of the mass contained in the inner regions of 
the halo. Lastly, all DM halos were constructed with no net 
angular momentum. 

The surface density profiles of the dwarf disks follow 
an exponential distribution in cylindrical radius R, while 
the vertical structure is mode led by consta nt-thickness, self- 
gravitating isothermal sheets ( Spit zerlfr942h 

p d (R,z) oc exp (-^) sech2 (j^) > ( 2 ) 

where Rj and Zd denote the radial scale length and the (sech 2 ) 
vertical scale height of the disk, respectively. The phase-space 
DF of the disk is fully determined once the disk velocity ellip- 
soid is specified. The radial velocity dispersion, <tr(R), is as- 
sumed to be exponential with (t\{R) = exp (— R / Rd) , where 
iTro denotes the central radial velocity dispersion. The disk az- 
imuthal dispersi on, ad>(R), is related to ctr (R) via the epicycle 
approximation (Bin ney & Trem aine 2008), while the vertical 
velocity dispersion, <r z (R), is set by the requirement that the 
adopted value of the scale height Zd is maintained in the total 
potential of the galaxy model. Here we parametrize the disk 
thickness as Zd/Rd and the disk mass as a given fraction, ma, 
of the halo mass, M/,. For simplicity, we also assume a con- 
stant value of <t S q = 10 kms" 1 for all our dwarf models. This 
choice results in a decreasing ratio of V Iot /a* with decreasing 
mass, where V rot and er* denote the stellar rotational veloc- 
ity and line-of-sight central velocit y dispersion, a trend that is 
indeed observed in dlrrs of the LG dMateoll 19981) . 

Our goal is to assess the degree to which the structure of 
late-type disky dwarfs can influence their tidal evolution in- 
side their host galaxies. For this reason, we generated a num- 
ber of dwarf galaxy models that differed in a number of impor- 
tant structural parameters, including the disk thickness, mass, 
and scale length, and the halo mass and concentration parame- 
ter. To investigate the impact of all these quantities on the tidal 
evolution of the dwarfs, we first constructed a "reference" or 
"default" dwarf galaxy model and subsequently initialized ad- 
ditional models by varying (increasing and decreasing) all rel- 
evant parameters in a systematic way, modifying only one at a 
time. Throughout the paper, we compare the effect of chang- 
ing a single parameter in three distinct dwarf models, namely 
the reference one and those with the largest and smallest val- 
ues of a given parameter. We also stress that we do not at- 
tempt any explicit rescaling of our dwarf galaxy models with 
redshift to account for the cosmic epoch at which the dwarfs 
wer e accreted by their hosts, as was done in previous stud- 
ies dMaver et al.l 1200 lag) . On the one hand, this allows us 
to avoid any uncertainties regarding the applicability of such 
scalings to dwarf galaxies. On the other hand, our system- 
atic parameter survey aims to address different choices for the 
basic structural parameters of the progenitors of present-day 
dSphs. As such, it should account to at least a certain extent 
for the wide diversity of dwarf galaxy structures expected at 
various cosmological epochs. 

In total, we constructed 1 1 high-resolution numerical mod- 
els of rotationally-supported dwarfs, which we denote Dl- 
Dll. The initial structural parameters of these models are 
listed in Table Q] (see § !2.4l for details on how the parameters 
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TABLE 1 

Initial Structural Parameters of the Dwarf Galaxy Models 







m d 






M h 


n 


Rd 




Vmax 


f max 


M/L 












Model 


zd/Rd 


(Mh) 


X 


c 


(1O 9 M ) 


(kpc) 


(kpc) 


2 


(kms- 1 ) 


(kpc) 


(M Q /L Q ) 


Vrot/f* 


P 


A 2 


b/a 


c/a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


Dl 


0.2 


0.02 


0.040 


20 


1 


1.29 


0.41 


3.93 


19.8 


2.07 


29.5 


2.93 


0.15 





1.00 


0.12 


D2 


0.1 


0.02 


0.040 


20 


1 


1.29 


0.41 


3.90 


19.9 


2.12 


30.4 


3.04 


0.29 





1.00 


0.06 


D3 


0.3 


0.02 


0.040 


20 


1 


1.29 


0.41 


3.91 


19.8 


2.19 


31.1 


2.85 


0.02 





1.00 


0.18 


D4 


0.2 


0.01 


0.040 


20 


1 


1.29 


0.41 


7.60 


19.4 


2.44 


64.3 


2.85 


0.21 





1.00 


0.12 


D5 


0.2 


0.04 


0.040 


20 


1 


1.29 


0.41 


2.05 


21.0 


1.69 


14.2 


3.06 


0.03 





1.00 


0.13 


D6 


0.2 


0.02 


0.024 


20 


1 


1.29 


0.25 


2.25 


19.8 


1.99 


27.2 


2.71 


0.05 





1.00 


0.11 


D7 


0.2 


0.02 


0.066 


20 


1 


1.29 


0.66 


6.32 


19.8 


2.32 


35.8 


3.06 


0.27 





1.00 


0.12 


D8 


0.2 


0.02 


0.040 


10 


1 


2.58 


0.41 


2.97 


16.5 


4.54 


42.4 


2.07 


0.27 





1.00 


0.10 


D9 


0.2 


0.02 


0.040 


40 


1 


0.65 


0.41 


5.00 


24.8 


1.19 


32.6 


3.78 


0.13 





1.00 


0.13 


D10 


0.2 


0.02 


0.040 


20 


0.2 


0.76 


0.24 


7.95 


11.6 


1.29 


31.2 


1.38 


0.58 





1.00 


0.12 


Dll 


0.2 


0.02 


0.040 


20 


5 


2.21 


0.70 


2.14 


32.0 


3.52 


25.9 


4.38 


-0.51 





1.00 


0.12 



NOTES. — The quantities in each column are as follows. Column 1: Dwarf galaxy model. Column 2: Scale height of the disk of the dwarf in units of the disk 
radial scale length, Rd- Column 3: Mass of the disk of the dwarf in units of the mass of the dwarf halo, M;,. Column 4: Spin parameter of the DM halo of 
the dwarf used to determine the scale length of the dwarf disk (see text for details). The corresponding disk scale lengths are listed in column 8. Column 5: 
Concentration parameter of the DM halo of the dwarf. Column 6: Mass of the DM halo of the dwarf in units of 10 9 Mq . This parameter is equivalent to the 
cosmologically-motivated virial mass. Column 7: Scale radius of the DM halo of the dwarf in kpc. Column 8: Radial scale length of the disk of the dwarf 
in kpc. Column 9: Toomre stability parameter of the dwarf disk at R = 2.5Rd- Column 10: Maximum circular velocity of the dwarf in kms -1 . Column 11: 
Radius at which the maximum circular velocity occurs in kpc. Column 12: Mass-to-light ratio of the dwarf in units of Mq/Lq. Column 13: Ratio of stellar 
rotational velocity to one-dimensional stellar velocity dispersion of the dwarf. Column 14: Anisotropy parameter of the stellar distribution of the dwarf. Column 
15: Amplitude of the m = 2 Fourier component of the surface density distribu tion of the dwarf stars. Column 16: Axis ratio b/a of the stellar component of the 
dwarf. Column 17: Axis ratio c/a of the stellar component of the dwarf. See jj |2.4| for details on how the parameters listed in columns 12-17 are calculated. Note 
that the entries in these columns are computed within r max . 



corresponding to columns (12)-(17) of this table are calcu- 
lated). Particular emphasis should be placed on column 4 of 
this table which lists the dimensionless spin parameters of the 
dwarf DM halos, A. This parameter is a measure of the to- 
tal energy content of a DM halo stored in rotation and it is 
defined as A = J\E\ '/ 2 /GM 5 / 2 dPeeblesll 19691) . where J is the 
total angular momentum, E is the binding energy, and M is 
the mass of the halo. Although our halos are non-rotating, 
we still employ the parameter A to assign scale lengths to 
the dwarf disks. This is because the disk scale length is not 
a free parameter in our m odeling, but rath er is derived via 
the semi-analytic model of Mo et al.l {1998) for the structure 
of disk galaxies in the ACDM paradigm 6 . According to this 
model, the baryons settle into a rotationally supported struc- 
ture whose scale length is determined by the mass, spin pa- 
rameter and concentration of the DM halo, and the fraction 
of mass and angular momentum in baryons relative to that of 
the halo. Assuming that the specific angular momentum of 
baryons is conserved during their infall and that the halo and 
baryons start with the same specific angular momentum, Rd is 
uniquely det ermined in our mode l s by M ; ; , A, c, and nid- We 
note that the lWidrow & D ubinski (2005) method for building 
disk galaxies does not explicitly take into account the adia- 
batic contraction of the ha lo in response to the slo w accumula- 
tion of the baryons (e.g., Blumenth al et al.l [l986). Therefore, 
in order to be consistent, we derive the values for the disk 
scale lengths without considering this effect. In summary, the 



6 In reality, the IMo et al] (1998) formalism and its assumptions for the 
formation of galactic disks may be inappropriate at the scales of dwarf galax- 
ies. This is due to the greater importance of thermal over rotational sup- 
port in such low -mass systems which modify the simple analytic scalings of 
IMo et all fT998l) . Nonetheless, given that our purpose is not to derive ex- 
act scale lengths for our dwarf galaxies but simply to determine the degree 
to which the tidal evolution of d isky dwarfs is affected by the sizes of their 
disks, following Mo et al. 1 1998) is reasonable and does not bias our results 
in any way. 



values of A that we discuss throughout this study do not reflect 
the angular momentum content of the DM halos, but rather 
serve the practical purpose of enabli ng us to derive th e values 
of Rd in the dwarf disks according to lMo et ail d!998l) . 

Our reference dwarf model Dl is characterized by the fol- 
lowing values for the adopted set of parameters: Zd /Rd = 0.2, 
m d = 0.02, A = 0.04, c = 20, and M h = 10 9 M Q . The details 
of how we varied these values in our simulation campaign 
are described in each relevant subsection of § [4] In what fol- 
lows, we motivate the choices for the default values of the var- 
ious parameters in dwarf model Dl. Our choices are mainly 
guided by studies of the properties of cosmological halos. 

Using a constraine d simulation of the LG, 
iKlimentowski et all d2010 ) studied the distribution of subhalo 
masses at the time of infall onto the primary as a function 
of redshift. These authors found that a significant fraction 
of satellites that were accreted by their hosts since z < 2 
and survived until the present time had mass es of « 10 9 M Q 
(see Figure 8 of IKlimentowski et all [2010). Assuming a 
concordance ACDM model, the median concentration value 
for a z = cosmological halo at this mass sc ale is c sa 20 
(e.g.. iBullock et aLllioofbt iMaccio et alj|2007n . In addition, 
both observational and theoretical evidence suggests that 
dwarf galaxies are not formed as thi n disks, but rather are 
born as thick, puffy systems (e.g., Dalc anton et alj l2004t 
iKaufmann etal]l2007l) . This is a consequence of the greater 
importance of feedback processes and turbulent motions in 
dwarf galaxies. We take this expectation into account by 
conservatively adopting Zd/Rd = 0.2, instead of the typical 
value Zd/Rd = 0.1 employed throughout the literature which 
would be more ap propriate for massive disk galaxies (e.g., 
iKregel et al.ll2002f) . 

As has been established by a large number of stud- 
ies, the distribution of halo spin parameters in TV-body 
simulations is well described by a log-normal distribu- 
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tion, with median values of A me d ~ 0.04 and dispersions 
of ax ~ 0.5 (e.g. . iBullock et alJ ^OOTi LShaw et al.l [2005 
IMaccio et alj|2007t iBett et alJl2007tlMacci6 et alj|2008l) Ac- 
cording to [M accio et al. (2007) (see also lBett et al.ll2007l and 
IMaccio et al.ll2008l) . the distribution of halo spins shows no 
dependence on halo mass. We stress that this conclusion 
does hold for halos with masses of the order of - 1O 9 M 
that we consider here (A. Maccio 2010, private communi- 
cation). Lastly, we choose the default value for the disk 
mass fraction of our dwarfs equal to mj = 0.02. This value 
is much lower than the universal baryon fraction but quite 
typical for pr esent-day low surface brightness (LSB) or din- 
galax ies (e.g. JJimenez et al.ll2003l:lGeha et al.ll2006tlOh et all 
2008), and is also in agreement with results of h ydrodynam- 
ical s i mulations of dw arf galaxy formation (e.g., Tass is et all 
120031: iGovernato et alJUoiOh . 

For each dwarf galaxy model, we generated an A^-body re- 
alization containing a total of 2.2 million particles (N/, = 10 6 
DM particles and Nd = 1 .2 x 10 6 disk particles). The gravita- 
tional softening was set to e/, = 60 pc and e ( / = 15 pc for the 
particles in the two components, respectively. The process 
of transforming a disky dwarf into a dSph via tidal stirring is 
fairly compl ex and depends on a number of subtle dynamical 
effects (e.g jMaver et al.l2001al) . It is thus important to estab- 
lish both the quality and the sufficient resolution of our dwarf 
galaxy models. For this reason, we evolved all dwarf galax- 
ies in isolation for a period of 10 Gyr. These test simulations 
revealed that the dwarfs retained their equilibrium configura- 
tion within the adopted force resolution over the timescales of 
the experiments. Therefore, our models should be largely un- 
affected by both two-body relaxation and artificial numerical 
heating of the dwarf disk through interactions with the mas- 
sive particles of the dwarf halo. The same simulations also 
confirmed the stability of all dwarf galaxy models against bar 
formation. Thus, any significant bar growth identified in the 
dwarfs during the course of the numerical experiments should 
be the result of the tidal field of the host galaxy, rather than a 
consequence of noise present in the initial conditions. All of 
these precautions imply that our dwarf galaxy models should 
indeed be adequate to resolve the generic tidal evolution of 
late-type disky dwarfs and to elucidate their transformation 
via tidally-induced dynamical instabilities and impulsive tidal 
heating. 

2.2. Primary Galaxy Model 

Another improvement we introduce in the present study 
is the fact that we employ self-gravitating primary galaxies 
as opposed to static host potentials adopted in the majority 
of earlier related investigations. The motivation behind this 
choice is twofold. First and most importantly, live primary 
galaxies can trigger instabilities in the disks of the dwarfs, 
that otherwise may not develop, and which could influence the 
dynamical a nd morphological evolu tion of the dwarf galaxies 
themselves (Weinberg & Blitz 2006). Second, by represent- 
ing the host galaxies as a distribution of interacting particles, 
we enable the dwarf galaxies to suffer dynamical friction. As 
a result of a progressively decaying orbit, the mass loss and 
tidal stripping experienced by the dwarfs will be enhanced 
with obvious consequences for their dynamical evolution in- 
side the host. 

For simplicity, we assume a single primary galaxy with the 
present-day structural pro perties of the MW In partic ular, we 
employ model MWb of Wid row & DubinskH d2005l) . which 
satisfies a broad range of observational constraints for the 



MW galaxy. Specifically, the exponential stellar disk has a 
mass of M D = 3.53 x 10 10 M Q , a radial scale length of R D = 
2.82 kpc, and a s ech 2 scale height of zd = 400 pc. The bulge 
which follows the Hernquist ( 1990) density profile has a mass 
and a scale radius of M B = 1.18 x 1O 1O M and as = 0.88 kpc, 
respectively. The DM halo has an NFW profile with a tidal ra- 
dius of R H = 244.5 kpc, a mass of M H = 7.35 x 10 n M o , and 
a scale radius of r# = 8.82 kpc. 

The simulations reported here use No = 10 6 particles in the 
disk, N B = 5 x 10 5 in the bulge, and N H = 2x 10 6 in the DM 
halo of the host galaxy, and employ a gravitational softening 
of tu = 50 pc, es = 50 pc, and = 2 kpc, respectively. The 
choice for the fairly large softening in the DM particles of the 
primary galaxy was motivated by our desire to minimize dis- 
creteness noise in the host potential. Such noise may lead to 
spurious two-body heating between the excessively massive 
halo particles of the primary galaxy and those of the dwarf 
disk, potentially interfering with the interpretation of our re- 
sults. In order to confirm the adequacy of our choice for 
we placed the reference dwarf model Dl at rest at various dis- 
tances inside the host galaxy and monitored the changes of 
its basic properties (e.g., surface density, velocity dispersions, 
thickness) as a function of time. The evoluti on o f the struc- 
tural parameters within radii of interest (see § !2.4l below) was 
found to be fairly small, specifically < 20%, during a period 
of several Gyr. This indicates that the chosen numerical pa- 
rameters are indeed appropriate to suppress the effect of two- 
body heating between the halo particles of the primary galaxy 
and those of the dwarf disk. 

2.3. Description of Tidal Stirring Simulations 

Our simulation campaign comprised 17 numerical experi- 
ments, which we denote R1-R17, of the tidal interactions be- 
tween rotationally-supported dwarf galaxies and their hosts. 
Table [2] provides a summary of all simulations we performed 
in this study. The various dwarf galaxy models are placed on 
bound orbits inside the primary galaxy. Except for varying 
the initial structure of the disky dwarf galaxies, we also inves- 
tigated the degree to which their tidal evolution is affected by 
the orbital parameters and, in particular, on the sizes and the 
eccentricities of the dwarf orbits. 

The size of the orbit is expected to be crucial for the ef- 
ficiency of tidal stirring, since the tidal force exerted by the 
primary galaxy depends strongly on the distance from its cen- 
ter. To explore this effect, we conducted a set of simulations 
where the default dwarf galaxy model Dl is placed on three 
orbits with different sizes. Although the orbits of LG dwarfs 
are currently poorly constrained observationally, their current 
distances, which give an indication of the apocenters of their 
orbits, coupled with studies of the orbital properties of cosmo- 
logical halos, can be used to inform our choices. Indeed, us- 
ing a c osmological simulation of the LG, Klimentows ki et all 
(2010) studied the orbital distribution of present-day satellites 
that were identified inside the virial radius of their hosts and 
had completed at least one orbit around them. These au- 
thors found that a significant fraction of this satellite pop- 
ulation had apocentric distances between r apo ps 250 and rj 
85 kpc, with typical values o f r apo w 125 kpc (see Figure 6 
of iKlim entows kiet alJ l2010l) . Adopting R vlI » 250 kpc as 
the v irial radius of the MW halo at z = (e.g., Kly pin et alJ 
l2002h . the previous numbers would roughly correspond to 
^vir> ^vir/3, and Ryjr/2. The value r apo = 125 kpc also ap- 
proxi mately matches the v irial radius of the MW halo at z ~ 1 
(e.g.. IWechsler et aT]|2002l) . A critical reader may note that 
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TABLE 2 
Summary of Simulations 



Simulation Dwarf Model (kpc) 
(1) (2) (3) 



(kpc) r apo /r pcri (deg) z d /Rd 
(4) (5) (6) (7) 



m d M h 
(M h ) X c (1O 9 M ) 
(8) (9) (10) (11) 



Rl 


Dl 


[2S 


2s 


5 


45 


0.2 


02 


040 


20 


R2 


Dl 


85 


17 


5 


45 


0.2 


0.02 


0.040 


20 


R3 


Dl 


250 


50 


5 


45 


0.2 


0.02 


0.040 


20 


R4 


Dl 


125 


12.5 


10 


45 


0.2 


0.02 


0.040 


20 


R5 


Dl 


125 


50 


2.5 


45 


0.2 


0.02 


0.040 


20 


R6 


Dl 


125 


25 


5 





0.2 


0.02 


0.040 


20 


R7 


Dl 


125 


25 


5 


90 


0.2 


0.02 


0.040 


20 


R8 


D2 


125 


25 


5 


45 


0.1 


0.02 


0.040 


20 


R9 


D3 


125 


25 


5 


45 


0.3 


0.02 


0.040 


20 


RIO 


D4 


125 


25 


5 


45 


0.2 


0.01 


0.040 


20 


Rll 


D5 


125 


25 


5 


45 


0.2 


0.04 


0.040 


20 


R12 


D6 


125 


25 


5 


45 


0.2 


0.02 


0.024 


20 


R13 


D7 


125 


25 


5 


45 


0.2 


0.02 


0.066 


20 


R14 


D8 


125 


25 


5 


45 


0.2 


0.02 


0.040 


10 


R15 


D9 


125 


25 


5 


45 


0.2 


0.02 


0.040 


40 


R16 


D10 


125 


25 


5 


45 


0.2 


0.02 


0.040 


20 


R17 


Dll 


125 


25 


5 


45 


0.2 


0.02 


0.040 


20 



0.2 
5 



NOTES . — Columns 3-5 refer to the initial orbital parameters of the disky dwarfs. Columns 7-1 1 list the initial structural parameters of the dwarf galaxy models. 
The quantities in each column are as follows. Column 1: Abbreviation for the tidal stirring simulations. Column 2: Dwarf galaxy model. Column 3: Apocentric 
distance of the orbit of the dwarf in kpc. Column 4: Pericentric distance of the orbit of the dwarf in kpc. Column 5: Eccentricity of the orbit of the dwarf. 
Column 6: Inclination of the disk of the dwarf with respect to the orbital plane in degrees. Column 7: Scale height of the disk of the dwarf in units of the disk 
radial scale length, Rj. Column 8: Mass of the disk of the dwarf in units of the mass of the dwarf halo, M;,. Column 9: Spin parameter of the DM halo of the 
dwarf used to determine the scale length of the dwarf disk (see text for details). The corresponding disk scale lengths are listed in column 8 of Table[T] Column 
10: Concentration parameter of the DM halo of the dwarf. Column 1 1 : Mass of the DM halo of the dwarf in units of 10'A/q . This parameter is equivalent to the 
cosmologically-motivated virial mass. 



the characteristic orbits of the z = subhalo population that 
we have considered may not be representative of those at the 
time of satellite infall onto the primary. This is especially true 
for systems that are accreted at high redshift (z > 1). Given 
the complexity of halo formation in a cosmological context, 
this is a valid concern. While addressing this issue is cer- 
tainly beyond the scope of the present paper, we note that at 
least dynamical friction should not have a major effect, if any, 
in altering the orbits of satellites wit h masses of the or der of 
~ 10 9 M Q that we have adopted here (IColpi et al.ll 19991) . 

Motivated by the previous discussion, we adopt r apo = 125, 
85, and 250 kpc as the apocentric distances of the dwarfs 
in the experiments with different sizes of orbits. We refer 
to these simulations a s "R l", "R2", and "R3", respectively, 
and discuss them in § 13.11 Fixing the eccentricity of the or- 
bits will determine the corresponding pericentric distances. 
Specifically, we adopt an eccentricity of r apo /V per i = 5, close to 
the median ratio of apocentric to perice ntric radii found in 
cosmological TY-bod y simulations (e.g., Ghigna et alj |1998t 
Diem and et al .1120071) . Interestingly, the pericentric distance 
of the tightest orbit with r per i = 17 kpc roughly corr esponds to 
that in ferred for the Sagittarius dwarf galaxy (e.g.. lLaw et al.1 
120051) . In addition, the pericentric distance of the most ex- 
tended orbit with r per i = 50 kpc is similar to that of the Large 
Magellanic Cloud (LMC), provided that the L MC has just 
crosse d the pericenter of its orbit around the MW (Bes la et ail 
120071) . Given the typical parameters associated with exper- 
iment Rl, throughout this work we refer to it as the "refer- 
ence" or "default" simulation that we use as the basis for the 
comparison with other experiments. 

Apart from the size of the orbit, we also investigate the 
impact of the orbital eccentricity on the tidal transformation 



of the dwarfs. Cosm ological simulations (e.g..|Ghigna et al.l 
Il998t iDiemand et al.ll2007t iKlimentowski et al.ll2010l) . indi- 
cate that satellite orbits range from nearly circular (r apo /V pel j ss 
1) to highly eccentric (r ap0 /r per i > 10). On high-eccentricity 
orbits, the effective duration of the tidal shock becomes so 
short that the response of the system is prone to be impul- 
sive rather than adiabatic and, as a result, ti dal heating is par- 
ticularly efficient (e.g. JGnedin et al.1ll999h . To this end, we 
considered two additional simulations in which we placed the 
dwarf galaxy model Dl on a highly and a mildly eccentric 
orbit, respectively. In particular, we considered orbits with 
eccentricities that were by a factor of 2 larger (r apo /V pel -i = 10) 
and smaller (r ap0 /r per i = 2.5) compared to that of Rl . We refer 
to these experi men ts as "R4" and "R5" (see Table |2} and dis- 
cuss them in § 13.21 In order to meaningfully compare the re- 
sults of this set of simulations, we kept the apocentric distance 
of the reference experiment Rl constant (r apo = 125 kpc). The 
desired value of eccentricity is then achieved by simply vary- 
ing the pericentric distance, r pel -j. This modeling results in a 
pericentric distance of r per i =12.5 kpc in simulation R4, which 
is the smallest in all experiments we performed in the present 
study. 

Initial conditions for the tidal stirring experiments were 
generated by building models of dwarf galaxies and placing 
them at the apocenters of their orbits. In the coordinate sys- 
tem chosen to describe the simulations, the orbital plane is the 
xy plane, and the center of mass of the combined system of 
dwarf and primary coincides with the coordinate origin. We 
stress that we did not impose any truncation in the density dis- 
tribution of the dwarfs (i.e.,, at the corresponding Jacobi tidal 
radius) to reflect their placement in the host tidal field. In ad- 
dition, the dwarf galaxies are not grown adiabatically in their 
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orbits, but rather are introduced in the simulations directly. 
We have explicitly checked that this set of assumptions does 
not bias our results. Indeed, the initial starting positions of the 
dwarfs correspond to a rather low density in the host galaxy. 
This results in a fairly weak tidal perturbation that does not 
affect the very inner regions of the dwarf, which constitute 
the primary target of the present study. 

According to t h e LG cosmological simulation of 
iKlimentowski et aTl (1201 0l) . most satellites that survived 
until the present time and whose masses at the epoch of infall 
onto the primary were in the range that we consider here 
(around 10 9 M Q ), entered their hosts at z < 2. Motivated 
by this finding, we follow the tidal evolution of the dwarfs 
inside their host galaxies for 10 Gyr. In the great majority 
of our experiments, the alignments of the internal angular 
momentum of the dwarf, that of the primary disk and the 
orbital angular momentum were all mildly prograde and 
equal to 45°. We discuss the implications of this choice in 
§ 15.61 but emphasize at the outset that our results should not 
be affected by any strong coupling of angular momenta. 

Lastly, all numerical experiments discussed in this work 
were performed w ith PKDGRA V, a multistepping, parallel, 
tree A^-body code ([Stadel 2001). In all experiments, we set 
the base time step to be equal to 1% of the dynamical time at 
the half-mass radius of the dwarf model and allowed the indi- 
vidual particle time steps to be at most a factor of 2 20 smaller. 
The time integration was performed with high enough accu- 
racy such that the total energy was conserved to better than 
0.5% in all cases, which is adequate for the type of study 
that we undertake in this paper. As we have already stated, 
our main goal is to investigate the dynamical evolution of the 
dwarfs which is driven by the response of their inner regions 
to the tidal shocks. The total energy contained in the inner 
parts of our dwarf models is a few percent of that of the entire 
system, so the accuracy of the energy conservation must be at 
least comparable to that in order to resolve meaningfully the 
dynamics of the region of interest. 

2.4. Parameters of Relevance 

In order to illustrate the evolution of the intrinsic, global 
properties of the dwarf galaxies as they orbit inside their 
hosts, we calculated a set of parameters as a function of time. 
Specifically, we examined the degree to which the mass, size, 
shape, and kinematics of the dwarfs is affected subject to the 
strong tidal field of the primary galaxy. 

For each simulation output, we first constructed the circular 
velocity profile of the dwarf galaxy, V c (r) = [GM(< r)/r\ l l 2 , 
where M(< r) is the spherically-averaged total mass profile 
about its center. As a measure of how the mass of the dwarf 
is affected at any stage, we employed the maximum circular 
velocity, V max , which reflects the mass distribution in the in- 
ner regions and is a commonly adopted quantity throughout 
the literature. The attractive feature of V max is that it is well 
defined and allows us to overcome d ifficulties regardin g the 
determination of the tidal radius (e.g. jReade t al. 2006b|). 

In order to demonstrate how the size of the dwarf galaxy is 
affected by tides, we investigated the ev olution of the radius 
at whic h V max occurs, r max . Following Klimentows ki et all 
(2009a), we chose r max as the characteristic scale where we 
computed the remaining properties of the dwarfs. This radius 
is large enough to include a significant fraction of the mass of 
the dwarf at all times, while at the same t ime is small enough 
not to be a ffected by the tidal tails (see IKlimentows ki et ail 
120071 l2009bl) . Due to the fact that tidal stripping truncates 



the dwarf galaxies at increasingly smaller physical radii, this 
characteristic radius decreases with time, from a few kpc ini- 
tially down to a fraction of a kpc at the end of the evolution. 

Because of the copious amounts of DM in dSphs, another 
quantity of particular interest is the mass-to-light ratio, M/L. 
Since the bulk of the stella r populations in these galaxies is old 
(t > 8- 10 Gyr; see, e.g., Grebel 2000), we chose to calculate 
the M/L ratio by conservatively adopting a fixed stellar mass - 
to-light ratio of (M/L)* = 3M /L (e.g. JSchulz eTaT]|2002h . 
We note that while such values for (M/L)* should be appro- 
priate for the present time, they may not be necessarily repre- 
sentative for the whole evolution, as the dwarfs might undergo 
periodic bursts of star formati on triggered by tid al compres- 
sion at pericentric passages (May er et alj|2001al) . However, 
for dwarf galaxies accreted by their hosts as early as 10 Gyr 
ago, as assumed here, most of the gas should be quickly 
stripped by ram pressure aided by the cosmic ionizing back- 
ground radiation (M ayer et al.ll2.007h . As a result, star forma- 
tion should cease soon after the dwarfs are accreted by the 
host galaxy. In addition, using the stel l ar population synthe- 
sis models of iBruzual & Chariot! dl993l) . lMaver et alJd2001al) 
showed that even when a burst of star formation occurs at first 
pericentric approach, the final (M/L)* ratio would correspond 
to that of an old stellar population after 7-8 Gyr of evolution 
inside the primary potential. This is because the effect of the 
newly-formed population on the (M/L)* ratio and on the color 
of the dwarf galaxy diminishes quite rapidly. The above dis- 
cussion suggests that the assumption of (M/L)* = 3M Q /L Q in 
the present study is fairly reasonable. 

Obviously, one of the most representative properties of 
dSphs is their spheroidal shape. We quantified the shape of 
the stellar component by calculating the moments of the iner- 
tia tensor for all stars within r max and deriving principal axis 
ratios b/a and c/a, where a, b, and c denote the major, in- 
termediate, and minor axis of the stellar distribution, respec- 
tively. Another interesting quantity which can illuminate the 
dynamical and morphological evolution of the dwarf galaxies 
inside their hosts, is the amplitude of the m = 2 Fourier com- 
ponent of the stellar surface density distribution, A2, given by 



1 

N 



2 ><t>i 



(3) 



Here ; denotes the two-dimensional cylindrical polar angle 
coordinate of particle j projected onto the xy plane (along the 
shortest axis) and the summation is performed over all stars 
within r max . The motivation behind studying the evolution 
of A2 is twofold. First, it can be used as an auxiliary mea- 
sure of the shape of the stellar com ponent since the tr iaxiality 
parameter, T = (a 2 — b 2 ) / (a 2 — c 2 ) (Fra nx et"aT1ll991l) . can be 
very noisy for systems close to spherical symmetry. We note, 
however, that A2 can still be ambiguous as its value is close 
to zero for both nearly spherical and disky shapes. Second, 
as we discuss below, in the majority of cases we studied, the 
strong tidal shocks at pericentric passages trigger the forma- 
tion of bars in the disks of the dwarfs. The parameter A2 will 
allow us to determine when a bar forms and to quantify the 
strength of the bar instability. 

For the purposes of the present study, we designate the for- 
mation of a bar when the amplitude of the m = 2 Fourier com- 
ponent satisfies A2 > 0.2 between two consecutive pericen- 
tric passages, namely for one full orbital period of the dwarf 
around its host. These conditions ensure that any tidally- 
induced bars that we detect in our simulations will be rela- 
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tively strong (see, e.g.,|Debattista et al. 2006 for even weaker 
bars in isolated disk galaxies). They also guarantee that any 
temporary increase of A%, when the dwarfs are strongly elon- 
gated by tidal forces at the pericenters of the orbit, will not be 
ascribed to a true bar instability. 

dSphs are low-angular momentum systems, as highlighted 
by their low ratio of rotational velocity to line-of-s ight central 
velocity dispersion, y ro t/c* < 1 (e.g jMateoll 1998b . It is thus 
important to investigate the evolution of the kinematics in the 
tidally-stirred disky dwarfs. To this end, after determining the 
directions of the principal axes of the stellar component in 
each output, we introduced a spherical coordinate system (r, 
9, <fi) so that the z-axis is along the shortest axis of the stellar 
distribution and the angle cf> is computed on the xy plane. We 
measured the kinematics of the dwarfs within r max by calcu- 
lating the rotational velocity around the shortest axis V lot = Va, 
and the dispersions o>, cq and cr^ around the mean values. 
We combined these dispersions into the one-dimensional stel- 
lar velocity dispersion parameter, cr* = [(of + crg + cr?)/3] 1 / 2 , 
which we adopted as a measure of the amount of random 
motions in the stars. We note that throughout this paper 
we compare the derived values of er* with those of observed 
dSphs. Defining cr* as an average over three directions makes 
such comparisons more meaningful, since the observed one- 
dimensional stellar velocity dispersions depend on the ran- 
dom line-of-sight. 

Another illuminating quantity related to stellar kinemat- 
ics is the anisotropy parameter, j3 = l-<7 2 /2er 2 , where o> 
and a, = (o'g + o'V) 1 ^ 2 denote the radial and tangential veloc- 
ity dispersion, respectively. This quantity describes the de- 
gree of anisotropy of the velocity distribution by indicating 
the amount of radial (fi = 1) versus circular (fi — > — oo) or- 
bits in the stellar population. We stress that this parameter 
cannot be measured directly from obse rvations, but only de- 
termined by dynamical modeling (e.g., lLokasI 120021) . Given 
that a well established range for the velocity anisotropy of ob- 
served dSphs does not yet exist, our simulations can provide 
useful predictions for the possible values of j3 in such systems. 

2.5. Classification Criteria 

The goal of the present study is to determine the conditions 
under which systems with the properties of dSphs can be pro- 
duced via the tidal interactions between disky dwarfs and their 
host galaxies. Observed dSphs possess a unique set of kine- 
matic and structural properties that we can utilize in order to 
perform comparisons with simulated dwarfs. In particular, we 
imposed two criteria for establishing the formation of a dSph 
in our simulations. 

The first criterion is related to the shape of the stellar com- 
ponent. Specifically, only simulated dwarf galaxies whose fi- 
nal states are characterized by c/a > 0.5 may be regarded as 
dSphs. The motivation behind this choice stems from the fact 
that most classic dSphs (e.g., Fornax, Draco, Leo I, Tucana, 
Sextans, Carina) have projec ted ellipticities, e = l-b/a, in 
the ra nge 0.1 < e < 0.5 (e.g.. lMateolfl998HMcGaugh & Wold 
2010), or equivalently projected axis ratios b/a > 0.5, where 
b and a denote the minor and major axis of the stellar dis- 
tribution, respectively. Given that an elongated shape can- 
not appear more elongated in projection, the condition that 
intrinsic c/a > 0.5 in the simulated dwarfs (n.b. by defini- 
tion b/a > c/a) ensures that the projected b/a will also sat- 
isfy b/a> 0.5 in any possible projection, in accord with the 
values inferred for observed dSphs. We note that the intrin- 



sic axis ratios of the dwarf remnants do not vary significantly 
with radius. Thus, they are appropriate for comparisons with 
observed ellipticities that are not computed within r max , but 
rather estimated eith er at the outer parts of the stellar distri- 
bution (iMateoll 19981) or within some observationally-defined 
radius (e.g., half-light radius). 

The second classification criterion is related to the kine- 
matics of the stellar component. In particular, we classify as 
dSphs only those systems which, after 10 Gyr of evolution 
inside the primary galaxy, exhibit V mt /a* < 1 in accordance 
with the res ults of observational studies of classic dSphs (e.g., 
Mateo 1 19981) . We stress that the values we compute for V rot 
and ct* are fairly close to what an observer would measure us- 
ing line-of-sight velocities, making the comparison with ob- 
servations reasonable. Indeed, at r max the rotation is close to 
the maximum value and can be readily measured if the line- 
of-sight is perpendicular to the rotation axis. Moreover, the 
velocity dispersion profiles of our remnants do not strongly 
vary with radius (as well as with line-of-sight) suggesting that 
our derived values for <r* should not differ significantly from 
those of the central velocity dispersion commonly used by ob- 
servers. As a result of all these facts, the condition Vrot/o* < 1 
that we impose here guarantees that V rot /a^ < 1 will also be 
satisfied for a significant fraction, if not all, of random lines- 
of-sight. 

Overall, for the purposes of the present study, the above 
classification scheme is adequate to establish the formation of 
a dSph in our simulations. A more direct comparison between 
theoretical and observational measurements for the properties 
of dSphs will be performed in a companion paper (Lokas et 
al. 2010, in preparation). 

3. EFFICIENCY OF TIDAL STIRRING AND ORBITAL PARAMETERS 
OF THE PROGENITOR DISKY DWARFS 

In this section, we gauge the dependence of the efficiency 
of tidal stirring on the orbital parameters of the progenitor 
rotationally-supported dwarfs, and in particular, on the size 
and eccentricity of their orbits. 

3.1. Size of the Orbit 

We first investigate the degree to which the tidal evolution 
of a disky dwarf is affected by the size of its orbit. The size 
of the orbit should be crucial for the outcome of tidal stir- 
ring since the tidal force depends strongly on the distance of 
the dwarf from the center of the host galaxy. In order to as- 
certain this, we considered two additional simulations (R2, 
R3) in which we placed the dwarf galaxy model Dl on orbits 
with sizes that were different from that of the reference ex- 
periment Rl. Specifically, in this set of simulations, we kept 
the eccentricity constant, but varied the default apocentric and 
pericentric distance by the same factor, leaving all other sim- 
ulation parameters unchanged (see Table |2j. Figure Q] shows 
the orbital trajectories of the dwarf galaxies from the respec- 
tive center of their host as a function of time for experiments 
Rl, R2, and R3. In Figure 12 we prese nt th e time evolution 
of the relevant parameters discussed in §|23](V max , r max , b/a, 
c/a, Vrot/c*, P, M/L, and A2) for the same simulations. 

In order to highlight some general trends in the evolution 
of these quantities with time, we first focus on the reference 
experiment Rl. As a result of the continuous action of the 
host galaxy tidal field, the overall structure of the orbiting 
dwarf is altered as it experiences a gradual decrease of its 
mass and physical size. This is reflected in the evolution of 
'"max and V max , respectively, which are both reduced systemat- 
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FIG. 1 . — Distances of the orbiting dwarf galaxies from the centers of their hosts as a function of time. Results are presented for simulations R1-R5, where the 
dwarfs are characterized by different orbital parameters. 



ically as a function of time. The preferential stripping of the 
mass from the outer, low-density regions of the dwarf galaxy 
leads to the adjustment of the circular velocity profile in such 
a way that both r max and V max decrease. The temporary but 
substantial increase of r max at the first pericentric approach 
can be attributed to the fact that the dwarf is being deformed 
by strong tidal forces for the first time and, as a result, its 
density distribution becomes very extended. This increase in 
'"max is not observed at subsequent pericentric passages as the 
dwarf galaxy becomes progressively more compact during its 
evolution and thus the ability of the tidal field to distort it is 
reduced. As we show below, in cases where the dwarfs are ini- 
tially more concentrated (as a result of either being embedded 
in halos of higher concentration or hosting more compact or 
massive disks) or are on wider orbits and experience weaker 
tidal forces, this increase in r max does not occur. 

The evolution of V max shows two interesting features. First, 
Vmax decreases significantly near the pericenters of the or- 
bit, where the intensity and variation of the time-dependent 
tidal force is the strongest, and the tidal shocks occur. This 
drop in V max is associated with mass loss from all regions of 
the dwarf galaxy. Because tidal shocks act on a very short 
timescale, even the inner regions of the dwarf, which are char- 
acterized by the highest densities and thus the shortest dy- 
namical times, respond impulsively to the external tidal per- 
turbation and mass is tidally stripped directly from regions 
within r max . However, V max also decreases because a fraction 
of the particles that are removed from larger radii are, in fact, 
on eccentric orbits with large apocentric distances. Although 
these particles spend most of their time outside of r max , they 
still contribute to the total mass in the inner regions. These 
are preferentially DM particles, since stars are confined to the 
central region of the dwarf by construction. After experienc- 
ing the tidal shock, it takes a few hundred million years, i.e., a 
multiple of the dynamical time of the system, for the dwarf to 
readjust to a new equilibrium state with a new value of V max . 

Second, V max remains remarkably constant between peri- 
centric approaches as the dwarf responds adiabatically to the 
weak intensity and variation of the tidal force. Overall, tides 
act nearly impulsively at pericentric passages and have lit- 
tle influence on the dwarf galaxy between pericentric pas- 
sages. This general behavior of V max and r max is observed in 
all our simulations, an d is in agreement with results reported 
in earlier studies (e.g.. | Hayashi et a l. 2003; Kazantzid is et al.l 
l2004bHKravtsov et al.ll2004tlPenarrubia et alJl2008l) . We note 



that, in the reference experiment Rl, the dwarf looses ~ 90% 
of its initial mass within r max during its orbital evolution but 
still survives as a bound entity. For such substantial mass loss, 
Vmax and r max decreased by a factor of ~ 1.7 and ~ 4, respec- 
tively. This finding has interesting implications as it suggests 
that tidal stripping affects the evolution of the r max - V max re- 
lation expected in the ACDM cosmological model (see also 
lKravtsovll201(ih . 

Focusing on the evolution of the M/L ratio, which we also 
estimate within r max , Figure [2] shows that it decreases mono- 
tonically as a function of time. This is a consequence of two 
facts. First, the stellar distributions of the dwarfs are natu- 
rally less extended than their DM halos. This implies that, 
within r max , the particles with the smallest binding energies 
and longest dynamical timescales, namely the particles that 
are most susceptible to be unbound by the tidal shocks, belong 
predominantly to the extended DM halos. Second, as dis- 
cussed above, particles that are stripped from the outer parts 
of the dwarf but still contribute to the total mass in the inner 
regions belong preferentially to the DM component. For these 
reasons, DM is stripped more readily than the stars within 
'"max and this decreases the M/L ratio. In addition to these 
two processes, during the orbital evolution of the dwarfs, r max 
decreases mono tonic ally and moves towards the inner, stellar- 
dominated regions, resulting in a natural increase of the frac- 
tion of mass in stars within r max . 

The evolution of the Fourier component A2 shows a sud- 
den increase after the first pericentric approach which signi- 
fies the formation of a tidally-induced bar. The onset of the 
bar instability coincides with a drop of V^t/V*, whose value 
continuously decreases thereafter. Loss of angular momen- 
tum caused by the bar and simultaneous increase of the stellar 
velocity dispersion due to tidal heating lead to a final value of 
Vrot/c* ~ 0.55. As a result of the tidal shocks, the strength 
of the bar diminishes as a function of time and the dwarf be- 
comes progressively more spherical. After 10 Gyr of evolu- 
tion inside the host, the initially-disky stellar distribution is 
transformed into a spheroid with axis ratios of c/a rj 0.7 and 
b/aTnQ.9. 

In what follows, we shall concentrate on comparing the 
tidal response of the dwarf in the reference simulation Rl with 
that in experiments R2 and R3. Figure|2]demonstrates that the 
dependence of the evolution of r max , V max , and M/L ratio on 
the size of the orbit is dramatic. As expected, the tighter the 
orbit is, the stronger and more rapid the decrease in all param- 
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FIG. 2. — Comparison of the evolution of various parameters as afunction of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the sizes of their orbits. Results are presented for the default simulation Rl and for experiments R2 and R3. The description of the simulations is presented in 
Table|2] For a fixed eccentricity, r a p /r pcl i, rotationally-supported dwarfs on tighter orbits with smaller pericentric distances, r pcl j, exhibit stronger tidal evolution 
inside their host galaxies and the efficiency of their transformation into dSphs is enhanced considerably. 



eters. Specifically, the dwarf galaxy in simulation R2 with the 
smallest pericentric distance looses ~ 99% of its initial mass 
within r max during its orbital evolution. In this case, V max and 
?"max decreased by a factor of ~ 2.7 and ~ 1 1 .5, respectively. 

In addition, Figure [2] illustrates two interesting features in 
the behavior of the dwarfs in experiments Rl and R2. First, 
the monotonic decrease of the M/L ratio is reversed for the 
tightest orbit in simulation R2 during the late stages of the 
evolution, indicating that the stars begin to be stripped more 
effectively compared to the DM. Second, the systematic de- 
crease of Vrot/c* is reversed temporarily for simulation Rl 
between the second and third pericentric passage. We post- 
pone the investigation of both issues for §|5] 

In the case of experiments Rl and R2, the evolution of the 
anisotropy parameter j3 reflects the transition to more radial 
orbits after the first pericentric approach, as the dwarfs change 
shapes to prolate spheroids. The stellar orbits tend to become 
more isotropic with time, especially for simulation R2 where 
the velocity distribution of the dwarf at the end is completely 
isotropic (J) pa 0). The velocity anisotropy of the remnant in 



experiment Rl is mildly radial with /3 pa 0.3. 

Interestingly, the most strongly perturbed dwarf in simula- 
tion R2 is almost entirely pressure-supported (Vm/v* ^5 0.1) 
and its final shape is also spherically-symmetric (b/a pa c/a pa 
1). Furthermore, while the dwarf in experiment R2 develops 
a strong bar, its counterpart in simulation R3 does not. The 
least evolved dwarf in experiment R3 remains oblate for the 
whole time except for a short period after the first pericentric 
approach. According to the criteria described in § 12.51 the fi- 
nal dwarfs in simulations Rl and R2 would be classified as 
dSphs but the remnant in experiment R3 would not as it has a 
disky shape and it is still dominated by rotation (VJ-ot/o* > 1). 
It is also important to highlight the strong association between 
bar formation and the formation of dSphs reported in this set 
of simulations. Such a connection is observed in all but one 
of our experiments. 

The final structural parameters of the simulated dwarfs af- 
ter 10 Gyr of evolution inside the primary galaxy are listed in 
Table [3] Special emphasis should be placed on columns 10 
and 11 of this table. Column 10 refers to whether a tidally- 
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TABLE 3 
Summary of Results 





Vmax 


T max 


M/L 
















Torb 


Simulation 


(kms- 1 ) 


(kpc) 


(Mq/Lq) 


Vrot/CT» 





A 2 


b/a 


c/a 


Bar Formation 


Classification 


(Gyr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Rl 


11.9 


0.53 


10.7 


0.55 


0.34 


0.06 


0.90 


0.73 


yes 


dSph 


2.09 


R2 


7.3 


0.18 


19.2 


0.08 


0.00 


0.00 


1.00 


1.00 


yes 


dSph 


1.28 


R3 


17.3 


1.22 


16.9 


1.46 


0.52 


0.14 


0.82 


0.26 


no 


non-dSph 


5.40 


R4 


7.2 


0.24 


15.3 


0.03 


0.03 


0.03 


0.95 


0.94 


yes 


dSph 


1.81 


R5 


16.0 


0.98 


14.0 


1.25 


0.54 


0.06 


0.88 


0.36 


no 


non-dSph 


2.50 


R6 


12.2 


0.51 


10.3 


0.24 


0.51 


0.04 


0.95 


0.73 


yes 


dSph 


2.09 


R7 


12.4 


0.55 


10.7 


0.23 


0.48 


0.16 


0.80 


0.74 


yes 


dSph 


2.09 


R8 


12.7 


0.50 


10.3 


0.74 


0.42 


0.10 


0.86 


0.66 


yes 


dSph 


2.09 


R9 


12.0 


0.55 


11.8 


0.61 


0.23 


0.04 


0.94 


0.77 


yes 


dSph 


2.09 


RIO 


10.4 


0.55 


18.8 


0.56 


0.19 


0.04 


0.93 


0.80 


yes 


dSph 


2.09 


Rll 


16.0 


0.49 


6.7 


0.40 


0.48 


0.28 


0.67 


0.58 


yes 


dSph 


2.09 


R12 


14.6 


0.41 


6.7 


0.56 


0.45 


0.29 


0.67 


0.53 


yes 


dSph 


2.09 


R13 


10.4 


0.56 


22.3 


0.26 


0.13 


0.03 


0.95 


0.91 


yes 


dSph 


2.09 


R14 


7.3 


0.40 


6.5 


0.31 


0.15 


0.01 


0.98 


0.89 


yes 


dSph 


2.10 


R15 


20.3 


0.67 


23.3 


1.18 


0.49 


0.07 


0.88 


0.52 


no 


non-dSph 


2.08 


R16 


7.1 


0.36 


11.6 


0.62 


0.33 


0.04 


0.94 


0.73 


no 


dSph 


2.14 


R17 


19.0 


0.77 


9.2 


0.70 


0.23 


0.09 


0.87 


0.68 


yes 


dSph 


1.88 



NOTES . — Columns 2-9 refer to the final structural parameters of the simulated dwarfs. The quantities in each column are as follows. Column 1 : Abbreviation for 
the tidal stirring simulations. Column 2: Maximum circular velocity of the dwarf in kms~'. Column3: Radius at which the maximum circular velocity occurs in 
kpc. Column 4: Mass-to-light ratio of the dwarf in units of Mq /Lq . Column 5: Ratio of stellar rotational velocity to one-dimensional stellar velocity dispersion 
of the dwarf. Column 6: Anisotropy parameter of the stellar distribution of the dwarf. Column 7: Amplitude of the m = 2 Fourier component of the surface 
density distribution of the dwarf stars. Column 8: Axis ratio b/a of the stellar component of the dwarf. Column 9: Axis ratio c/a of the stellar comp onent of 
the dwarf. Column 10: Formation of a tidally-induced bar during the orbital evolu tion of the dwarf galaxies according to the criteria discussed in § 12.41 Column 
11: Classification of the final systems according to the criteria imposed in § 12.51 Column 12: Orbital time in Gyr defined as the average time elapsed between 
consecutive apocentric passages. Note that the entries in columns 4-9 are estimated within r maj[ . 



induced bar was formed in the stellar distribution of the dwarf 
galaxy during the course of its evolution inside the host. Col- 
umn 1 1 indicates whether a dSph was produced as a result of 
the tidal interaction between the progenitor disky dwarf and 
the primary galaxy. The entries in column 6 demonstrate that 
the values of f3 in the dSph remnants range from isotropic to 
mildly radial (0 < (3 < 0.5). Given that a well established 
range for the velocity anisotropy of observed dSphs does not 
yet exist, a direct comparison with observations is not possi- 
ble. However, we note that the magnitude of the radial veloc- 
ity anisotropy in our dSphs is rather similar to that observed 
in the outer p arts of CDM halos for med in cosmological sim- 
ulations (e.g. JCole & Lacevl lT996). Lastly, the orbital times 
listed in the last column of Table [3] are defined as the average 
time elapsed between consecutive apocentric passages. 

Figure [5] shows the surface density maps of the final stel- 
lar distributions of the dwarfs in simulations Rl, R2, and R3. 
Results are presented for projections onto the xz, yz, and xy 
planes and particles are color-coded on a logarithmic scale, 
with hues ranging from orange to white indicating increasing 
stellar density. This figure visually confirms the conclusions 
advanced above regarding the degree of morphological trans- 
formation experienced by the dwarfs in in this set of experi- 
ments. 

3.2. Eccentricity of the Orbit 

In this section, we explore the extent to which the tidal evo- 
lution of a rotationally-supported dwarf is influenced by the 
eccentricity of its orbit. As we discussed earlier, this effect is 
worth investigating as the effective duration of the tidal shock, 
and thus the response of the system to t he tidal perturbation , 
depends on the orbital eccentricity (e.g.. lGnedin et alj|1999l) . 



To this end, we conducted two additional simulations (R4, 
R5) in which we placed the dwarf galaxy model Dl on orbits 
with eccentricities that are larger (r apo /r per i = 10) and smaller 
(fapoAperi = 2.5) by a factor of 2 from that of the reference 
experiment Rl (see Table©. Simulations R4 and R5 should 
be viewed as corresponding to a highly and a mildly eccen- 
tric orbit, respectively. Figure Q] shows the orbital trajectories 
of the dwarf galaxies from the centers of their hosts in these 
experiments. 

Figure [4] shows the main results related to this set of simu- 
lations. As expected, this figure demonstrates that for a fixed 
apocentric distance, higher eccentricity orbits induce a much 
stronger transformation in the orbiting dwarfs. The evolution 
of the various parameters in this set of experiments proves 
to be analogous to the one presented in Figure [2] The dwarf 
galaxies in simulations R4 and R5 evolve similarly to those 
in R2 and R3, respectively, and the same strong link between 
bar formation and the formation of dSphs is reported. Indeed, 
the disky dwarf in experiment R4 develops a bar and the fi- 
nal product is a nearly spherical (b/a w c/a > 0.9), isotropic 
(J3 w 0.03), and non-rotating (V rot /er* « 0.03) dSph, with final 
properties akin to those of the remnant in simulation R2. The 
evolution of V max , r max , and mass within r max is also similar in 
the two cases. This is due to the fact that although the orbital 
time of the dwarf in experiment R4 is significantly larger, is 
has a smaller pericentric distance compared to that of R2. On 
the other hand, the dwarf in experiment R5 does not form a 
bar and remains oblate, radially anisotropic and rotating, anal- 
ogous to the one in simulation R3, despite the shorter orbital 
time and larger number of pericentric passages (see Figure[T)). 

A last interesting thing to note is that the bar instability that 
develops in the most strongly perturbed dwarf (R4) is weaker 
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FIG. 3. — Surface density maps of the final stellar distributions of the dwarfs in simulations Rl (upper panels), R2 (middle panels), and R3 (lower panels). 
In order to aid comparison, the first panel also includes the edge-on view of the initial disk. Particles are color-coded on a logarithmic scale, with hues ranging 
from orange to white indicating increasing stellar density. Local density is calculated using an SPH smoothing kernel of 32 neighbors. Results are presented for 
projections onto the xz (left columns), vz (middle columns), and xy planes (right columns), where the x, y, and z axes lie along the major, intermediate, and minor 
axes of the stellar distribution, respectively. Rl is our reference simulation. Experiment R2 with the shortest orbital time and the second smallest pericentric 
distance all our experiments (7" or b = 1.28 Gyr and r pcr j = 17 kpc) yields a spherically- symmetric (b/a RJ c/a ~ 1) and isotropic (fi 0) stellar system with 
negligible amounts of rotation (Vrot/f* J$ 0. 1) that would be classified as a dSph. The dwarf in simulation R3, whose orbit is characterized by the longest orbital 
time and largest pericentric distance (T or \, = 5.40 Gyr and r^, = 50 kpc), is not transformed into a dSph and remains disky even after 10 Gyr of tidal evolution 
inside the host galaxy. 



compared to that of the reference simulation, despite the fact 
that the pericentric distance is much smaller in the former 
case. The same behavior, although to a lesser degree, is ob- 
served in Figure |2] We speculate that this is a consequence 
of the very strong tidal forces in simulations R2 and R4. The 
tidal shocks may heat the newly formed bar so substantially 
that they actually cause the instability to rapidly dissipate. 

To summarize, the results of this and the previous section 
highlight that the effectiveness of the transformation into a 
dSph via tidal stirring depends crucially on the orbital param- 
eters of the progenitor disky dwarfs. The relative importance 



of the orbital time and the pericentric distance in this process 
will be addressed in §|5] 

4. EFFICIENCY OF TIDAL STIRRING AND STRUCTURAL 
PARAMETERS OF THE PROGENITOR DISKY DWARFS 

In this section, we explore how the initial structure of the 
progenitor disky dwarfs could influence the outcome of their 
tidal evolution inside the gravitational field of their host galax- 
ies. We begin by focusing on the parameters of the dwarf disk 
and discuss in turn the effect of th e disk inclination, thickness, 
mass, and scale length. In § 14.21 we turn our attention to the 
properties of the DM halo of the dwarf and investigate the ex- 
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FIG. 4. — Comparison of the evolution of various parameters as afunction of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the eccentricities of their orbits. Results are presented for the default simulation Rl and for experiments R4 and R5. The description of the simulations is 
presented in TablefS] For a fixed apocentric distance, r apo , rotationally-supported dwarfs on orbits with higher eccentricities, r a p /r per i, and thus smaller pericentric 
distances, r pcr ;, display stronger tidal evolution inside their host galaxies and the efficiency of their transformation into dSphs is increased substantially. 



tent to which the effectiveness of tidal stirring is affected by 
the halo concentration parameter and mass. 

4.1. Varying the Parameters of the Disk 
4. 1 . 1 . Disk Inclination 

We first investigate the extent to which the tidal evolution of 
a late-type disky dwarf is affected by the initial inclination of 
its disk with respect to the orbital plane. This is worth explor- 
ing as stronger alignments between the orbital angular mo- 
mentum of the dwarf and the internal angular mo mentum of 
its disk can result in more effective stripping (e.g jRead et al.l 
2006b). For this purpose, we performed two additional sim- 
ulations in which we placed the dwarf galaxy model Dl on 
the same orbit as in simulation Rl after changing the default 
inclination from i = 45° to i = 0° and to i = 90° (see Table©. 
Figure [5] contains the results pertaining to this set of simula- 
tions. 

Although the initial inclination of the dwarf disk is very dif- 
ferent, the size (in terms of r max ), the mass (in terms of Vm ax ), 



and the M/L ratio evolve fairly similarly. This is a conse- 
quence of two facts. First, DM halos were constructed with 
no net angular momentum. Thus, their stripping should be in- 
dependent of such considerations and proceed in exactly the 
same way. Second, the dwarfs are not stripped down to such 
small scales that the alignment between the stellar angular 
momentum and the orbital angular momenta wo uld have an 
import ant effect on the stripping of the stars (e.g., iRead et all 
2006b). Indeed, we find the decrease in mass within r max for 
both stars and DM to be independent of the initial disk incli- 
nation. 

It is important to note that Klime ntowski et all (l2009al) re- 
ported that disk inclination has a significant effect on the evo- 
lution of V max in similar tidal stirring experiments. Specif- 
ically, they found that V max decreased much more strongly 
in the i = 0° case compared to the i = 45° an d i = 90° cases 
(see Figure 10 of IKlimentowski et al.ll2009al) This was be- 
cause the DM halos in the IKlimentowski et al . (2009a) exper- 
iments were constructed with net angular momentum which 
was aligned with the angular momentum of the dwarf disk. 
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FIG. 5. — Comparison of the evolution of various parameters as afunction of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the inclinations of their disks with respect to the orbital plane, i. Results are presented for simulations Rl, R6, and R7. For the specific choices of i in this 
set of experiments, the tidal evolution of the rotationally-supported dwarfs inside their host galaxies and the efficiency of their transformation into dSphs depend 
very weakly on the disk inclination. 



Moreover, the stellar disks were much more extended com- 
pared to those of the present study and, as a result, the strip- 
ping of the dwarfs down to the smallest scales of stars did 
occur at some point during the evolutio n. In fact, the mass 
within r max of both stars and DM in the Klimentows ki et all 
(2009a) experiments decreased more in the i = 0° case lend- 
ing support to the previous arguments. 

The behavior of the bar strength amplitude A 2 in Figure [5] 
demonstrates that all dwarfs possess a tidally-induced bar af- 
ter the first pericentric passage. However, the bar appears to 
be much weaker in simulation R7, indicating that the tidal 
forces acting on the dwarf disk were smaller in this case. 
During subsequent pericentric passages the situation becomes 
progressively more complex as strong variations in the values 
of Ai are observed among the three simulations. At the end 
of the evolution, the tidally-induced bars are diminished in all 
cases. Investigating the reasons for the very different evolu- 
tion of the bars in this set of experiments is clearly beyond the 
scope of the present paper. However, it is interesting to note 
that recent targeted numerical experiments have highlighted 



the importance of stochasticity in the evolution of isolated 
disk galaxie s leading to macroscopic diffe rences in the evolu- 
tion of bars dSellwood & D ebattista 200 9j). Among the causes 
for this stochastic evolution. ISellwood & Debattista! d2009l) 
identified interference between multiple disk modes, ampli- 
fication of noise, bending modes, dynamical friction between 
the bar and the halo, and intrinsic chaos. Obviously, further 
analysis would be required to determine whether the observed 
behavior of the bars in this set of simulations can also be at- 
tributed to stochasticity and/or possibly to other considera- 
tions such as the different initial disk inclinations and reso- 
nances. 

The final systems in simulations R6 and R7 are character- 
ized by c/a pa 0.7 and V ro t/c* ~ 0.2 (see TableO and, there- 
fore, would be classified as dSphs according to the criteria 
described in § 12.51 While the evolution of kinematics and 
shape is, in broad terms, fairly similar among the three exper- 
iments, there are some interesting differences that are worth 
mentioning. For example, the rotation is lost more quickly in 
simulation R6 (z = 0°), where Vrot/c* drops almost to zero at 
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FIG. 6. — Comparison of the evolution of various parameters as afunction of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the thicknesses of their disks. Disk thicknesses are parametrized as Zd/Rd and results are presented for simulations Rl, R8, and R9. For the specific choices of 
Zd/Rd in this set of experiments, the tidal evolution of the rotationally-supported dwarfs inside their host galaxies and the efficiency of their transformation into 
dSphs are essentially independent of the disk thickness. 



~ 5 Gyr and remains constant thereafter. The substantial and 
rapid decrease of V ro t/0* in this experiment is explained by 
the presence of a very strong bar instability between ~ 1 and 
~ 5 Gyr. We also note that in experiment R7 the axis ratio b/a 
remains almost constant with time and significantly different 
from unity, so that the stellar component is triaxial rather than 
prolate at early times. 

4.1.2. Disk Thickness 

In this section, we investigate the dependence of the tidal 
evolution of a disky dwarf on the thickness of its disk. Such 
an investigation is important for two reasons. First, feed- 
back mechanisms and turbulent motions in dwarf galaxies 
should be effective in producing thicker systems. Second, the 
efficiency of tidal heating and the strength of bar instabili- 
ties that are both vital for tidal s tirring should be differe nt in 
thicker stellar distributions (e.g., Kazantzidis et al.ll2009h . To 
this end, we generated two additional dwarf galaxy models 
that were identical to Dl except for their scale height, which 
was chosen equal to Zd/Rd = 0.1 and Zd/Rd = 0.3. We re- 



fer to these models as "D2" and "D3", respectively (see Ta- 
ble [TJ. The choice for the larger thickness in model D2 re- 
flects the greater importance of pressure support in low-mass 
galaxies and was motivated by both observations as well as 
results of recent numerical simulations of the formation of 
isolated galaxies (Kauf mann et al.ll2007l) . including systems 
with y max s» 20kms _1 as our default dwarf model. On the 
other hand, model D3 with the thinner disk was constructed 
mainly for completeness. After building these two new dwarf 
models, we placed them on the same orbit as in simulation Rl 
and followed their tidal evolution inside the host galaxy. We 
refer to these experiments as "R8" and "R9" (see Tabled and 
present the relevant results in Figure [6] 

Overall, for the specific choices of initial thickness in this 
set of experiments, no significant difference in the evolution 
of the dwarfs can be discerned. Had we adopted thicker 
and/or thinner disks, such differences may have been more 
pronounced. The final properties of the systems in simula- 
tions R8 and R9 indicate that the remnants can be classified 
as dSphs (see Tabled. 
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Figure [6] demonstrates that r max , V max , and the M/L ratio 
evolve almost identically in the three experiments. A simi- 
larly strong bar is induced in all dwarfs after the first pericen- 
tric approach. While the thicker disk is initially more isotropic 
by construction, all three dwarfs exhibit nearly identical val- 
ues of j3 between the first and second pericentric approach. 
At later times, the evolution is somewhat different, but gen- 
erally the the stellar component of the thicker disk (R9) has a 
more spherical shape and is characterized by more isotropic 
orbits. All dwarfs end up mildly triaxial and very similar to 
the remnant in simulation Rl, but their paths to this state are 
slightly different. Occasionally, departures from the typical 
trend of decreasing rotational velocity occur (see Figure [TTJ 
but they happen at different pericentric passages for different 
experiments. We come back to this issue in §[3] 

4.1.3. Disk Mass 

In this section, we investigate the degree to which the tidal 
evolution of a rotationally-supported dwarf is affected by the 
mass of its disk. In order to ascertain this, we generated 
two additional dwarf galaxy models that were identical to Dl 
except for their disk mass, which differed by a factor of 2: 
rrid = 0.01 and mj = 0.04 7 . We refer to these models as "D4" 
and "D5", respectively (see Table [TJ. This range of mj val- 
ues is consistent with that inferred from exquisite measure- 
ments of the mass distribution of nearby dwarfs combining 
high resolution gas kinematics in the THINGS survey and 
deep Spitzer photometry (Oh et al. 2010, in preparation). 
It is also in agreement with results of hy drodynamical sim- 
ulations of dwarf ga laxy formation (e.g.. iTassis et al.l 12003b 
iGovernato et al.|[2010t) . After building these two new dwarf 
models, we placed them inside the host galaxy on the same 
orbit as in simulation Rl and followed their tidal evolution. 
We refer to these experiments as "RIO" and "Rll" (see Ta- 
ble[2]> and present their results in Figure|7] 

This figure demonstrates that the mass of their disk has a 
substantial effect on the tidal evolution of the disky dwarfs. 
Interestingly, although the initial values of r max are different 
and smaller for the more massive disks due to the more con- 
centrated mass distribution, at the end of the evolution they 
converge to nearly the same values. The opposite is true for 
the mass, as quantified by V max . Indeed, the least massive disk 
increases the susceptibility of the dwarf to tidal effects and the 
mass loss within r max is larger in this case. Overall, the results 
presented in Figure [7] indie ate that the masses of the dwarfs 
are differentiated as a result of the tidal evolution. 

The M/L ratios, which differed significantly in the begin- 
ning, evolve in a fairly different way. The most massive disk, 
which had the smallest value of M/L initially, is affected most 
weakly and its M/L ratio decreases only slightly. The weak 
evolution of the M/L ratio in this case can be attributed to the 
fact that the initial r raax is very small (and remains so until 
the end of the evolution). As a result, r max probes the central 
regions o f the dwarf, where the reasons that we have indi- 
cated in § 13. 1 I for the preferential stripping of DM over stars 
do not apply. Indeed, DM particles within this inner region 

7 We note that formally, according to the model of Mo et al. 11998), vary- 
ing the disk mass and keeping the other relevant parameters (M;,, A, and c) 
constant would result in different scale lengths for the dwarf disk. However, 
for consistency and in order to isolate the effect of a single parameter on the 
tidal evolution of the dwarfs, we adopt identical disk scale lengths in this set 
o f expe riments. The same convention applies to the simulations described in 
§ 14.2. ll in which we construct dwarf models with different halo concentration 
parameters. 



are very unlikely to be on orbits of high enough eccentricity 
to have apocenters larger than the tidal radius and be stripped. 
Moreover, the slopes of the density profiles of stars and DM 
are similar below ~ 1 kpc, implying that the average bind- 
ing energies of DM and stars do not differ significantly in this 
region. Interestingly, despite the different evolution, the hier- 
archy of M/L is preserved; the dwarf which has the highest 
M/L ratio initially also exhibits the highest M/L ratio in the 
end. 

The mass of the disk affects the evolution of the stellar 
shape substantially, and to some degree also influences that 
of the stellar kinematics. Owing to its higher self-gravity, 
the most massive disk (Rll) develops the strongest tidally- 
induced bar after the first pericentric approach, leading to a 
pronounced decrease in V Iot /a^. We note that the values of 
A2 remain the largest in this case until the end of the evolu- 
tion (f = 10 Gyr). As a result, the dwarf galaxy in simula- 
tion Rll retains a strongly prolate shape with final values of 
c/a 0.6 and b/a » 0.7. On the other hand, the dwarf with 
the least massive disk (R10) remains triaxial rather than pro- 
late for most of the time. Due to the significantly lower disk 
self-gravity in this case, the tidal shocks heat the disk much 
more effectively, producing a quite spherical stellar compo- 
nent in the end (c/a w 0.8, b/a~ 0.9). 

However, the interpretation of the evolution of kinematics 
is less obvious. While the least massive disk loses its rota- 
tion monotonically (see Figure fTTh. the evolution of its most 
massive counterpart is more chaotic with occasional strong 
increases of the rotational velocity, even after two consecu- 
tive pericentric passages (third and fourth). We note that this 
behavior is observed in only one more simulation (R12; see 
next section) but is not as strong as reported here. Overall, the 
values of c/a, in conjunction with the fact that Vrot/c* 0.6, 
indicate that the final systems in simulations R10 and Rll 
would be classified as dSphs (see TableO. 

4. 1 .4. Disk Scale Length 

This section addresses the extent to which the tidal evolu- 
tion of a disky dwarf is influenced by the radial scale length 
of its disk, Rd- For this purpose, we generated two additional 
dwarf galaxy models that were identical to Dl except for the 
scale lengths of their disks. We reiterate that the disk scale 
length is not a free parameter in our simulations, but rather 
is derived via t he semi-analytic galaxy formation model of 
iMo et all (119981) . Bee ause of this choice our method requires 
us to assign values to the halo spin parameter A in order to 
derive the scale length of the dwarf disk, despite the fact that 
the halos of our dwarfs are constructed with no net angular 
momentum. 

Recall that we chose A = 0.04 for our reference dwarf 
model Dl, close to the median value o f halo spins found in 
cosmologic al Af-b ody simulations (e.g . .lBullock et aLlfe OOla: 
IShaw etalJ 12005 iMaccio et al.1 120071: iBett et alj|2007l) . To 
investigate the effect of the disk scale length, we initialized 
two additional dwarf models with A = 0.024 and A = 0.066. 
These values correspond to the la deviations from the me- 
dian A, assuming a log-normal distribution with A me d = 0.04 
and <x\ =0.5, in agreement with the previous cosmological 
studies. We refer to these dwarf models as "D6" and "D7", 
respectively, and list their disk scale lengths in Table Q] In 
the corresponding experiments "R12" and "R13", we placed 
these dwarfs on the same orbit as in simulation Rl and fol- 
lowed their tidal evolution inside the host galaxy. Note that 
for consistency we have kept the thickness Zd /Rd in this set 




FIG. 7. — Comparison of the evolution of various parameters as a function of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the masses of their disks. Disk masses are parametrized as a given fraction, mj, of the halo mass, M;,, and results are presented for simulations Rl, RIO, and 
Rll. Rotationally-supported dwarfs with less massive disks exhibit stronger tidal evolution inside their host galaxies and the efficiency of their transformation 
into dSphs is augmented significantly. 



of experiments constant. As a result, the disk becomes larger 
(smaller) in both the radial and vertical direction when A is 
increased (decreased). The relevant results are presented in 
Figure [8] which demonstrates that the remnant dwarfs in this 
set of simulations would be classified as dSphs (see TableO. 

This figure shows that the scale length of the dwarf disk 
is a parameter that affects significantly the tidal evolution of 
the disky dwarfs. The initial value of r max is smaller for the 
more compact disk in simulation R12 and it remains so until 
the end of the evolution. This experiment does not show the 
characteristic temporary increase of r max at the first pericentric 
approach observed in the reference simulation Rl. This is be- 
cause the more compact disk diminishes the ability of the host 
tidal field to distort the inner regions of the dwarf galaxy. As 
expected, the mass loss within r max measured by the decreas- 
ing value of V max is more pronounced for the more extended 
disks with a larger value of A. This is simply a consequence of 
the fact that larger disk scale lengths make the potential well 
of the dwarf shallower and decrease the total binding energy 
in the inner parts. 



As expected, the dwarf with the most extended disk (R13) 
has a larger initial value of the M/L ratio within r max . Due 
to the dissimilar values of Rj, the stripping of the stars pro- 
ceeds differently in this set of experiments and, as expected, 
is much more effective in the case of the most extended disk. 
This differential stripping of the stars manifests itself in the 
evolution of the M/L ratio. Although the initial values of the 
M/L ratios are not very different, the final value of M/L is 
much larger in simulation R13. Interestingly, the M/L ratio 
in this experiment starts to increase during the intermediate 
stages of the evolution (f ~ 5 Gyr). Such a behavior is sim- 
ilar to that reported in § [3] for experiments R2 and R4, even 
though this transition occurred earlier there and the final M/L 
values obtained were slightly lower. We discuss this issue in 
§|5] Overall, Figure [8] indicates that the M/L ratios are differ- 
entiated as a result of the tidal evolution. 

Evidently, the size of the disk affects substantially the evo- 
lution of the stellar shape, and to some degree also influences 
that of the stellar kinematics. The behavior of the Fourier 
component A2 after the first pericentric passage is slightly 
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FIG. 8. — Comparison of the evolution of various parameters as afunction of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the radial scale lengths of their disks, Rj. Results are presented for simulations Rl, R12, and R13. The spin parameters of the dwarf DM halos, A, used to 
determine the disk scale lengths in this set of experiments are indicated in the labels (see text for details). The corresponding values of Rj for the three different 
choices of A are equal to Rj = (0.41,0.25,0.66) kpc, respectively (see Table[T}. Rotationally-supported dwarfs with more extended disks experience stronger 
tidal evolution inside their host galaxies and the efficiency of their transformation into dSphs is enhanced substantially. 



misleading. Indeed, both lower and larger A cases show sim- 
ilar values of this parameter, but for different reasons. While 
the most compact disk (R12) retains an oblate shape, the 
most extended one (R13) is triaxial, but neither shows such 
a strong bar mode as our reference experiment Rl. The sit- 
uation changes radically after the second pericentric passage. 
On the one hand, owing to its higher self-gravity, the low- 
A dwarf forms a fairly strong bar and remains prolate until 
the end of the evolution with final values of c/a « 0.5 and 
bjaK, 0.7. In fact, the remnant in this experiment exhibits the 
least spheroidal shape of all our classified dSphs, highlighting 
the resilience of the most compact disk to tidal heating. This 
system could be the analogu e of the most elongated dSph s 
in the LG such as Ursa Minor tllrwin & Hatzi dimitriou 19951) . 
On the other hand, tidal effects are quite efficient in the case of 
the high-A disk with the lowest self-gravity and cause its rapid 
transformation into a nearly isotropic (f3 « 0.1) and spherical 
stellar component (b / 'a « c / a> 0.9) that retains little rotation 

(VrotK < 0.3). 



4.2. Varying the Parameters of the Halo 
4.2.1. Halo Concentration 

In this section, we investigate the degree to which the tidal 
evolution of a disky dwarf is affected by the concentration pa- 
rameter of its DM halo, c. To this end, we constructed two 
additional dwarf galaxy models that were identical to Dl ex- 
cept for the halo concentration parameters. Recall that we 
chose c = 20 for our reference model Dl, close to the me- 
dian concentrati on value for z = cosmological halos of m ass 
~ 1O 9 M re.g. JBuUock etlDl2Wb1:lMacci6 et alJl2007h . To 
investigate the effect of the halo concentration on our results, 
we constructed two additional dwarf models with c = 10 and 
c = 40. These values roughly correspond to the 2er deviations 
from the median c = 20, in accordance with the previous cos- 
mological studies 8 . We refer to these dwarf models as "D8" 

8 We note that c = 40 is markedly incompatible with results of the distri- 
bution of DM in present-day LSB and dlrr galaxies through the modeling of 
rotation curves (see the recent review of de Blok 2010 for a discussion on this 
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FIG. 9. — Comparison of the evolution of various parameters as a function of time illustrating the dependence of the tidal transformation of disky dwarf galaxies 
on the concentration parameters of their DM halos, c. Results are presented for simulations Rl, R14, and R15. Rotationally-supported dwarfs embedded in less 
concentrated DM halos exhibit stronger tidal evolution inside their host galaxies and the efficiency of their transformation into dSphs is increased considerably. 



and "D9", respectively, and list their structural parameters in 
TableQ] In the corresponding experiments "R14" and "R15", 
we placed these dwarfs on the same orbit as in simulation Rl 
and followed their tidal evolution inside the host galaxy. Fig- 
ure [9] contains the results of these simulations. 

As expected, the less concentrated dwarf is characterized 
by a larger initial r max . As in the case of the experiments with 
different disk masses (R10, Rll), after the first pericentric 
approach the values of r max converge. Interestingly, at late 
times the hierarchy is even reversed (the dwarf with the most 
concentrated halo exhibits the largest r max ) but the differences 
are not significant. As demonstrated by the weakest decrease 
in the value of V max , the most concentrated halo increases the 
resilience of the dwarf galaxy to tides and mass loss. Since 
the higher concentration makes the potential well of the dwarf 
deeper, this effect is similar to that of the larger disk mass 
discussed in S I4.1.3I The effect on the M/L ratio is, however, 
the opposite; although the M/L ratio within r max was initially 

issue). Nevertheless, even though it is clearly inconsistent with observations, 
we decided to consider this limiting case simply for completeness. 



larger for the c = 10 case it soon becomes the smallest. This 
is because the DM particles are much more weakly bound in 
this case and, as a result, are stripped fairly efficiently. After 
a period of t = 10 Gyr inside the host galaxy, the M/L ratios 
of the dwarfs with the most and the least concentrated halo 
differ by a factor of ~ 4. The weak evolution of the M/L ratio 
in simulation R15 is of similar origin to that of the experiment 
with the most massive disk (Rl 1). 

The susceptibility to tides of the low-concentration dwarf 
in experiment R14 also manifests itself in the evolution of 
the shape and kinematics. The stellar component in this ex- 
periment evolves rapidly and becomes nearly spherical rather 
quickly with a final c/a w 0.9. The stronger effect of the tidal 
shocks in this case reflects the fact that the dwarf responds 
more impulsively to the tidal perturbation owing to its lower 
characteristic densities and, correspondingly, longer internal 
dynamical times. After the first pericentric passage, a strong 
bar develops in the dwarf disk (A2 > 0.3) and the Vtot/c* ratio 
starts to decrease, reaching a final value of V rot /<7* ~ 0.3. The 
structure of the remnant in this case would resemble that of 



20 



KAZANTZIDIS ET AL. 



observed dSphs. 

On the other hand, the stellar component of the high- 
concentration dwarf retains its disky shape, since it is less 
affected by tides. This is similar to the cases of simulations 
R3 and R5 where the dwarf galaxies experienced weaker tidal 
forces due to the larger pericentric distances. In experiment 
R15, the final axis ratio is higher (c/a « 0.5) but still the stel- 
lar component remains oblate. Most importantly, as demon- 
strated by the very low values of A-i throughout the evolu- 
tion, a bar does not form in this case. As a result, the an- 
gular momentum content of the dwarf does not decrease ef- 
ficiently and the Vrot/c* ratio remains larger than unity until 
the end of the simulation. The absence of the bar in this case 
is consistent with results of numerical experiments showing 
that high values of halo concentration suppress the forma- 
tion ofbarsJnjsoMeddisk galaxies embedded in CDM ha- 
los dMaver & Wadslevl2004l) . To summarize, experiment R15 
constitutes the only case where a rotationally-supported dwarf 
is placed on the reference orbit of simulation Rl and neither 
develops a bar nor yields a dSph. This lends further support 
to the fact that bar formation and the formation of dSphs are 
intimately linked in the context of the tidal stirring model. 

4.2.2. Halo Mass 

In this section, we explore the degree to which the tidal 
evolution of a disky dwarf depends on the mass of its DM 
halo, Mh- Such a dependence may be expected as various 
effects, including the strength of tidal shocks as well as dy- 
namical friction, will be affected by the halo mass. For this 
purpose, we generated two additional dwarf galaxy models 
that were identical to Dl except for the masses of their ha- 
los, which differed by a factor of 5: M h = 0.2 x 10 9 M Q and 
Mh = 5 x 10 9 M Q . This mass range roughly corresponds to 
the range of initial subhalo masses that were accreted since 
z < 2 and survived until the pres ent time in the LG simula- 
tion of IKlimentowski et al.l d2010h . We refer to these models 
as "D10" and "Dll", respectively (see Table [TJ. The moti- 
vation behind using the same values for nid, Zd/Rd, c, and 
A in these models as in the reference one (Dl) was imposed 
mainly by simplicity. However, we do note th at our choices 
are reasonable given th e scatter in r rid (e.g., iJimenez et alJ 
20031: iTassis et al.ll2003l) . z d /Rd (e.g., Dalc anton et al l 120041: 
Kaufmann et al.1 120071) . and c (e.g., Bull ock et al.l [2001b: 
Maccioet al.1 120071)" for galaxies in this mass range, and the 
fact that the distribut ion of halo spins shows no dependence 
on halo mass (e.g. jMaccio et al.ll2007l) . In the corresponding 
experiments "R16" and "R17", we placed dwarf models D10 
and Dll on the same orbit as in simulation Rl and followed 
their tidal evolution inside the host galaxy. Figure [TOlpresents 
the results of these simulations. 

As a result of our choices for the initial structural param- 
eters of the dwarf galaxies in this set of experiments, r max is 
larger for the dwarf with the most massive halo (R17). This is 
simply because r max m 2.16/?/,/c for the NFW profile, where 
Rh in our modeling corresponds to the cosmologically moti- 
vated virial radius of the system (/?/, w r v j r ). Figure [TOl shows 
that r max remains larger for the most massive dwarf until the 
end of the evolution. We emphasize that the pericentric pas- 
sages which are marked by sudden drops in V max do not occur 
at the same time for all three dwarfs, despite the fact that the 
orbital parameters were initially identical. Specifically, the 
orbit of the most massive dwarf becomes progressively much 
tighter with continuously reduced apocentric distances, lead- 
ing to a shorter orbital time (see Table |3J. This is due to the 



effect of dynamical friction, whose strength depends strongly 
on the mass of the moving body. 

Interestingly, after 10 Gyr of evolution inside the host, all 
three dwarf galaxies lose approximately the same fraction 
of their initial mass within r max (~ 90%) and their V max de- 
creases by approximately the same factor of ~ 1.7. This is 
intriguing given that mass loss is generically affected, among 
other things, by the depth of the potential well of the dwarf 
and the details of the dwarf orbit (e.g., the pericentric dis- 
tance) around the primary galaxy. We note that by construc- 
tion the three dwarfs have approximately equal average densi- 
ties within r max , and therefore comparable dynamical times at 
r mal . Using this fact coupled with some order-of-magnitude 
calculations it is possible to gain qualitative insight as to the 
relative evolution of V max and mass within r max is similar in 
simulations Rl, R16, and R17. 

According to the tidal approximation of the impulse (which 
should be valid in our simulations for the purposes of dimen- 
sional analysis), the energy injected at each pericentric pas- 
sage is given by AE oc M^ oA MR 2 V~2, where Mh os t is the mass 
of the host, M denotes the mass of the dwarf, R is a char- 
acteristic radius of the dwarf (related to r max ), and V re i is the 
relative v elocity between the two galaxies at the pe ricenter of 
the orbit (ISpitzerlll958tlBinnev & Tremaind l2008). Since, as 
noted above, the average density within r max is roughly con- 
stant among the dwarfs, we have R oc M 1 ' 3 in the three dif- 
ferent models. Given that Mh os t is the same and V re i is fairly 
similar in the three experiments considered here, it follows 
that AE oc M 5 / 3 . By virtue of the virial theorem, the energy 
content of the dwarf should scale as E oc GM 2 /R<xM 5 /\ and 
also E oc MV 2 , where V is a characteristic velocity of particles 
within the dwarf (related to y max ). Hence, the fractional in- 
crease in energy caused by the tidal shocks, AE /E, is roughly 
constant, explaining why the V max and M max of the dwarfs in 
simulations Rl, R16, and R17 evolve in a fairly similar way. 
We stress that similar conclusions would be reached had we 
used extensions of the standard impulse approximation appro- 
priate for perturbers with extended mass distributions as well 
as for satellites that move on either straight-path or eccentric 
orbits dGnedin et al.lll999h . This is because such extensions 
are characterized by analogous scalings of the relevant vari- 
ables in our calculations above. 

The M/L ratios are similar in all cases and decrease mono- 
tonically with time reaching a value of ~ 1OM Q /L in the 
end of the evolution. The behavior of kinematics shows an in- 
teresting feature. While the values of V r ot/c* and j3 were quite 
different for the three models initially, they become very sim- 
ilar after the first pericentris passage and remain so until the 
end of the evolution. The evolution of the shape shows a more 
clear trend with mass. After the first pericentric approach, the 
bar is stronger in the most massive dwarf compared to that of 
the default case and it remains so until the end of the evolu- 
tion. This simply reflects the higher self-gravity of the disk in 
simulation R17. On the other hand, while the least massive 
dwarf remains disky during the entire evolution, its axis ratio 
c/a increases significantly with a final value of c/a ps 0.75. 
The dwarf remnants in simulations R16 and R17 are charac- 
terized by c/a « 0.7 and V ro t/ c* < 0.7 (see Tabled and hence 
would be classified as dSphs. 

A last point of interest is that simulation R16 constitutes 
the only case where the formation of a dSph is not associ- 
ated with the presence of a tidally-induced bar. Indeed, the 
increase of the bar amplitude, A2, in this case is only tem- 
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FIG. 10. — Comparison of the evolution of various parameters as a function of time illustrating the dependence of the tidal transformation of disky dwarf 
galaxies on the masses of their DM halos, M/,. Results are presented for simulations Rl, R16, and R17. Owing to the interplay between the effects of dynamical 
friction and the depth of the potential well, rotationally-supported dwarfs embedded in DM halos of substantially different mass display similar tidal evolution 
inside their host galaxies and the efficiency of their transformation into dSphs is of comparable magnitude. 



porary and occurs at the first pericentric passage when the 
dwarf is strongly elongated by tidal forces. The absence of a 
bar is consistent with the fact that the disk of this dwarf ex- 
hibits the highest value of the Toomre stability parameter (see 
Table Q]). Bar formation causes rotational motion to be trans- 
formed into radial motion and thus, in conjunction with tidal 
heating which also increases the velocity dispersion, is inti- 
mately linked to the decrease of V rot / cr* r equired for the trans- 
formation to a pressu re-supported dSph dMaver et al.l l2001a; 
iDebattista et al.ll2~006l) . In the case of the least massive dwarf, 
the initial Vrot/c* within r max is rather low (Vrot/c* < 1 -5) and 
therefore the bar stage is not necessary to decrease it consid- 
erably to levels that are appropriate for dSphs (V m t/u* < 1). 
The decrease of V rot / er* in the dwarf of experiment R 1 6 occurs 
almost exclusively at pericentric approaches and is driven by 
the decrease in V rot (see FigurefTTI). The decrease in rotational 
velocity within r max is simply caused by tidal heating which 
perturbs the circular orbits of the disk stars converting part 
of the rotation (ordered motions) in the disk to random mo- 
tions. It is worth noting that in the reference simulation Rl, 



for example, the combination of tidal heating and the tidally- 
induced bar lead to a decrease in the V ro t/c* of the dwarf by 
approximately a factor of 5, illustrating the importance of bar 
formation in reducing the initial high values of V mt /a* appro- 
priate for disky dwarfs to those characteristic of dSphs. 

5. DISCUSSION 

We begin this section by investigating the V msix -<7* relation 
in our simulations, which is of particular interest for the miss- 
ing satellites problem. In § 15.21 and § 15.31 we explore the rea- 
sons for the increasing rotational velocity and the increasing 
M/L ratio, respectively, observed in some of our experiments. 
In § 15.41 we elucidate the fundamental orbital parameter that 
influences the transformation of disky dwarfs into dS phs and 
thus the efficiency of the tidal stirring mechanism. In § 15.51 we 
discuss some additional implications related to the findings of 
this work. Lastly, in § 15.61 we conclude with a few words 
of caution and a discussion of promising directions for future 
work that may lead to more conclusive statements about the 
role of environmental processes in shaping the nature of dwarf 



22 



KAZANTZIDIS ET AL. 





15 


o 




CD 




W 


1 








5 









15 






QJ 




(/] 

\ 


10 


E 






c: 
U 









15 









10 


tn 








5 









15 


O 






1 u 


CO 


















15 


o 




OJ 




t/} 

\ 


10 


t 






5 












1 1 1 1 1 1 1 


M M M ' - 
R12 : 


: , 1 


,1,1, 


1 , 1 , 1 , 1 , : 




2 3 4 5 6 7 
time (Gyr) 



R13 : 



9 10 



-V ,: rotational velocity 
■ <J + : ID velocity dispersion 



- 1 1 1 1 1 1 ' 


III- 

R5 : 






: , 1 , 1 , 1 , 


l , i , l , l , l , : 




- 1 1 1 1 1 1 ' 


i | i | ' | ' | ' | ' . 

R8 : 


M 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , - 





' I 1 I 


1 1 I 1 I 


' 1 1 - 








Rll : 









, i , : 



-'I'M 


M ' 1 ' 1 ' 1 ' 1 ' 1 M 
R14 : 


: M , 1 M 


»T 1 M , 1 M , 1 , 1 r= 



1 5 
10 
5 




15 
10 
5 





- M ' 


MM 


1 ' 1 ' 1 


1 ' 1 M 
R4 : 


1 , 1 


,1,1 


j . . Vt 


-i — m— ; — 3 




-MM 


1 


1 ' 1 ' 1 M 






R10 : 








: , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , - 







- M 


' 1 M ' 1 


1 ' 1 ' 1 ' 1 ' - 






R16 - 


-.1,1,1,1,1, ,1,1,1,- 



23456789 10 
Lime (Gyr) 



1 



23456789 10 
Lime (Gyr) 



FIG. 1 1 . — Evolution of the stellar rotational velocity, V ro t, (solid lines) and the one-dimensional stellar velocity dispersion, cr*, (dotted lines) of the dwarfs 
in simulations R1-R17. Parameters Vrot and cr* are computed within r max - The repeated action of tidal forces exerted by the host galaxy decreases the ratio of 
Vrot/f*, gradually transforming the initially rotationally-supported disky dwarfs to pressure-supported stellar systems dominated by random motions. 



galaxies in environments such as that of the LG. 



5.1. V m 



- cr* Relation 



Here we investigate the relation between the one- 
dimensional, central stellar velocity dispersion, cr*, and the 
maximum circular velocity, V max , in our simulated dwarf 
galaxies. The knowledge of this relation is crucial, par- 
ticularly in the context of the miss ing satellite problem 
dMoore et al .11 19991; iKlvpin et alll 19991) . as it reflects the link 
between observable properties and dark halo properties in 
dwarf galaxies. One of the main difficulties intrinsic to this 
problem lies in the fact that the mapping between Vmax and 
a*, the latter being a directly observable quantity, depends 
sensitively on the assumed density structure of the DM halo 
of the dwarf. For example, adopting an isothermal halo model 
with a flat circular velocity profile to estimate V max would 
yield V max ~ \f2o* (e.g.. iMoore et al.|[l999l) . while assum- 



ing an isotropic stellar velocity dispersion tensor would result 
in Vmax ~ v^cr* (e.g.. IKlvpin et al.l[l999l) . Given t he uncer- 
taintie s regarding the conversion of cr* to Vm ax (see lKravtsovl 
1201 Ol for a thorough discussion on this issue), it is interest- 
ing to check the general validity of these assumptions directly 
in our simulations. The relevant analysis is presented in Fig- 
ure[[2] 

The left panel of this figure shows the relation between 
the one-dimensional stellar velocity dispersion and the maxi- 
mum circular velocity V max in the final stages of all simulated 
dwarfs (f = 10 Gyr). We remind the reader that the veloc- 
ity dispersion profiles of our remnants do not strongly vary 
with radius and thus our dispersion measurements are close 
to the central values commonly used by observers. For the 
purposes of this presentation, we have also replaced cr* used 
throughout the paper by a'^ = [(u^ + al + a'^)/3] x l 2 , where 



On the Efficiency of the Tidal Stirring Mechanism for the Origin of dSphs 



23 




5 io 15 5 10 15 5 10 15 

o-'^km/s) ff' # (km/s) 



FIG. 12. — Relation between the maximum circular velocity, Vmax, and the one-dimensional, central stellar velocity dispersion, a',, of the dwarfs in our 
simulations. The latter quantity corresponds to the line-of-sight central velocity dispersion measured in the observations and includes the contribution of rotation 
(see text for details). Left: The data points indicate the values calculated in the final stages of the dwarfs (t = 10 Gyr) in experiments R1-R17, and the dotted line 
is intended to identify a slope of \/3 and is not a fit to the data. The formal best-fit slope differs by less than ~ 4% from \/3. Circles show results for specific 
experiments Rl, R2, R3, and R5 discussed in the right panel. The relation Vmax = %/3o"», adopted by Klvoin et al. 1 1999) to map the observed central stellar 
velocity dispersion of dSphs to the maximum circular velocities of their halos, is well justified within the context of the tidal stirring model. Right: Four examples 
of the evolution of individual dwarf galaxies on the V mlix -cr' t plane (solid lines). Results are shown for representative simulations Rl, R2, R3, and R5, where the 
same dwarf model (Dl) was placed on different obits. Rl is our reference experiment. R2 produces one of the most strongly evolved dwarfs in our simulations, 
while R3 and R5 correspond to the least tidally-transformed systems in our experiments. In all panels, the filled circles mark the end of the evolution inside the 
host (/ = 10 Gyr). As in the left panel, dotted lines correspond to slopes of V3. The dwarf galaxies exhibit fairly different types of trajectories depending on their 
orbital parameters and can move substantially on the V max — o"^ plane over the course of their evolution inside their hosts. 
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■■ 00 + V r o t . This is because rotation was found to be im- 
portant in some of the remnants, specifically those of simula- 
tions R3, R5, and R15 which do not produce dSphs, and thus 
it had to be included in the calculation of their velocity disper- 
sion. We stress that this modification is not necessary in the 
majority of cases where dSphs with minimal intrinsic rotation 
are formed. 

The panel demonstrates a strong correlation between 
and V max in all dwarf remnants regardless of their morphol- 
ogy and kinematics. This is not entirely unexpected as both 
of these quantities are measures of the depth of the potential of 
the dwarfs. Interestingly, the relation V max = V3a r , formally 
valid for a tracer stellar population with an isotropic velocity 
dispersion tensor, reprod uces the results of the sim ulations re- 
markably well (see also Kli mentowski et aLll2009al) . Indeed, 
the formal best-fit slope of the V max - er^ relation in our sim- 
ulations differs by less than ~ 4% from V3- Given the wide 
range of orbital and structural parameters that we adopted for 
the progenitor disky dwarfs in our experiments, this finding is 
particularly noteworthy. We stress that the relation between 
er, and V max is derived by directly probing the stellar kinemat- 
ics of the simulated dwarfs. It may thus be reliably used to 
map the observed central stellar velocity dispersions to halo 
maximum circular velocities in dwarf galaxies when address- 
ing the missing satellites problem. 

Overall, we conclude that moderate conversion factors 
between a* and V max such as those origin ally adopted by 
iMoore et ail (1 1999b and lKlvpinet all d 19991) to formulate the 
missing satellites problem were reasonable. Recently, using 
controlled simulations of subhalo evolution, Penarrubi a et alj 



(2008) have argued in favor of much larger conversi on factors 
in ob served MW satellites (ss 2-3; see, however. iKravtsovl 
2010 f or a discussion of possible caveats in the approach of 
Penarrubia et al. 2008). We stress that our results should at 
least hold for dwarf galaxies that have moderate M/L ratios 
as our final systems (e.g., Leo I or Fornax). Dwarfs that are 
much more DM dominated (e.g., Draco or Ursa Minor) may, 
in principle, exhibit a different relation between cr» and V max 
owing to their ( possibly) differen t formation and evolutionary 
histories (e.g., Mayer et al. 2 0071) . 

Furthermore, the relation between V max and er* reported 
here is valid for dwarfs that were originally embedded in 
cuspy NFW-like DM halos, like the ones employed in the 
present study, and would likely be different with other halo 
density structures. Indeed, if the halo density profiles in the 
progenitor disky dwarfs were shallower than NFW in the in- 
ner re gions that are probed by the stars (e.g.. lGovernato et al.1 
2010), the resulting circular velocity profiles will still be 
slowly rising within the luminous radius, and the conver- 
sion factor betwe en a* and Vm ax ca n be substantially larger 
(> 2-4: see. e.g.. lStoehr eTaIll2002l) . 

Our findings are corroborated by a number of studies. For 
example, applying Jeans modeling and varying the an i sotropy 
in the velocity distribution of stars, Kazantzidis et al. ( 2004b) 
found that the stellar kinematics of Fornax (er* ~ 1 1 - 
13kms _I ) and Draco (cr* ~ 8-lOkms" 1 ), can be repro- 
duced in halos with V max ^20-30 kms" 1 , which would im- 
ply cr* /Vmax ~ 0.3 — 0.65. Similar results were obtained by 
IZentner & Bullockl (120031) and t his range of V max was also 
confirmed by lMadau et al.l (120081) using the cosmological Via 
Lactea simulation of the formation of a MW-sized halo. An- 
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FIG. 13. — Surface density map of the stellar distribution of the dwarf 
in simulation Rll. Particles are color-coded on a logarithmic scale, with 
hues ranging from yellow to blue indicating increasing stellar density. Local 
density is calculated using an SPH smoothing kernel of 32 neighbors. The 
distribution of stars is projected onto the orbital plane and is visualized at the 
third pericentric approach (t =5.25 Gyr from the start of the simulation). The 
solid line indicates the direction towards the center of the host galaxy, located 
at (x,y) = (0,0). A bar is evident in the stellar component and the rotational 
velocity of the stars, V ro t, is anticlockwise in this image. The tidal forces 
point in the same direction as the rotational velocity vectors of the stars in the 
dwarf, resulting in an increase of V ro t- 



other investigation of Draco by lLokas etail (120051) . which 
considered the contamination of the stars by tidal tails and 
included the fourth velocity moment in the analysis, found 
er„ ~ 8 kms" 1 and V max ~ 16 kms " 1 , aga in in agreement with 
our results. Lastly, Walker et al. (2009) have recently ana- 
lyzed the observed velocity dispersion profiles of a number of 
classic dSphs. Treating the anisotropy of the stellar velocity 
distribution as a free parameter, these authors estimated the 
y max of the host halos in the range ~ 10-25 kms" 1 , similar to 
those that would be obtained assuming the relation between 
Vmax and er» reported here. 

The right panel of Figure[T2]illustrates four examples of the 
evolution of individual dwarfs on the Vmax — plane. Results 
are presented for representative experiments Rl, R2, R3, and 
R5, where the same dwarf galaxy (Dl) was placed on differ- 
ent obits, but the general trends are confirmed in all simula- 
tions. The combination of tidal stripping and loss of angular 
momentum dominating over tidal heating dictates the charac- 
teristic evolution toward lower values of V max and a'^ seen in 
all cases. The panel also shows that the dwarf galaxies ex- 
hibit different trajectories depending on their orbital parame- 
ters and can move substantially on the Vmax _ plane during 
the course of their evolution inside their hosts. Both facts im- 
ply that the final V ma x _ &' t relation in our simulations is not a 
consequence of the initial conditions. 

5.2. Increasing Rotational Velocity 

In the majority of cases we studied, the rotational veloc- 
ity of the stars within r ma x, V rot , decreased systematically with 
time leading from the initial rotationally-supported dwarfs to 
pressure-supported stellar systems dominated by random mo- 
tions. There are instances in the evolution, however, where 
Vrot increases at pericentric passages. As seen in Figure [TT] 
this occurs in simulations Rl, R9, and R17 (after the sec- 



ond pericentric approach), R8 (after the third pericentric ap- 
proach), and Rl 1 and R12 (after the third and fourth pericen- 
tric passage). 

Investigating for the reason behind the temporary increases 
of the rotational velocity we checked the orientation of the 
elongated stellar component at the pericentric approaches at 
which such increases occur. We found that the typical orienta- 
tion of the bar is similar to that illustrated in Figure [13] where 
we see the stellar component of simulation Rl 1 at the third 
pericentric approach projected onto the orbital plane. The 
rotational velocity of the stars is anticlockwise in this plot. 
Therefore, the tidal forces which act along the direction to- 
wards the host galaxy (solid line) point in the same direction 
as the rotational velocity vectors of the stars in the dwarf, ef- 
fectively speeding up the rotation. Note that the case shown 
in Figure[l3]is the only instance in the entire evolution of the 
whole suite of 17 experiments where the total angular mo- 
mentum of the stars within r ma x increases significantly. In all 
other cases of increasing V rot the angular momentum does not 
increase because it is controlled not only by the rotational ve- 
locities of the stars but also by r max , which decreases at each 
pericentric passage. 

Lastly, we note that the observed increases in V 10t are al- 
ways associated with changes in the shape of the stellar com- 
ponent from a more prolate to a more triaxial configuration, 
or equivalently with an increase of the difference between the 
axis ratios b/a and c/a. This is most evident in experiment 
Rll, where the prolate shape with b/a w c/a « 0.4 at the 
third pericentric passage transforms into a triaxial one (with 
b/a > c/a). A similar phenomenon occurs at the subsequent 
(fourth) pericentric approach. These changes are also accom- 
panied by changes in other parameters. As verified in Figure|2] 
and Figures l3TIT0l when the rotational velocity (or V rot /(7*) in- 
creases, there is always a drop in the anisotropy parameter f3 
and the bar strength amplitude A2. This means that at these in- 
stances the shape of the stellar component becomes less pro- 
late and the orbital structure of the stellar distribution changes 
as the stellar orbits become less radial. 

5.3. Increasing M / L Ratio 

In most of our simulations, the M/L ratio computed within 
'"max decreases monotonically as a function of time for the rea- 
sons we discussed in § 13.11 In three of our 17 experiments, 
however, the M/L ratio starts to increase at some point during 
the evolution. This happens most prominently in simulations 
R2 and R4, which are characterized by the smallest pericentric 
distances and produce the most strongly tidally stirred dwarfs. 
Although much smaller, such an increase in the M/L ratio 
is also observed in experiment R13 where the stellar disk is 
more extended due to the larger initial value of the radial scale 
length, Rd, adopted in this case. We stress at the outset that 
these increases of the M/L ratio are rather moderate and thus 
cannot account for the exceptional DM content in some of the 
classic dSphs such as Draco and Ursa Minor. This implies that 
the very high values of M/L inferred for some dSph galaxies 
may not be a consequence of tidal evolution, but could be the 
result of either the formation process of their progenitors or 
other mechanisms that affected the baryonic mass fraction of 
dSphs. For example, the presence and subsequent removal 
of gas by the combination of tides and ram pressure has been 
sho wn to be vital in producing systems with extreme DM con- 
tent dMayer et alJl2006ll2007l) . 

Figure[l4]investigates the reason behind the increasing M/L 
ratio. This figure shows the relative mass loss, \AM\/M, of 
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stars and DM separately in simulations R2, R4, and R13. For 
comparison, we also present the same quantities in our refer- 
ence experiment Rl. To avoid noise in the data, the values of 
\AM\/M were assigned at pericenters by taking the masses 
within r max at the following and preceding apocenter, namely 
\AM\/M = \M M -Mj\/Mj, where i is the apocenter number. 
Figure [14] demonstrates that in all four simulations the rela- 
tive mass loss of the DM component is initially much larger 
than that of the stars, and that the former always dominates 
in experiment Rl. In simulations R2, R4, and R13, however, 
the stars begin to be stripped more effectively at some point 
during the late stages of the evolution, leading to the increase 
of the M/L ratios reported in Figures [2] @1 and [8] 

The reason for this behavior can be attributed to the fact that 
in these three cases the dwarf galaxies are stripped down to the 
scales where the alignment between the angular momenta of 
the stars and the orbital angular momenta of the dwarfs start to 
have an important effect on stellar stripping. In experiments 
R2 and R4, with the smallest pericentric distances, this is due 
to the very strong tidal forces. In simulation R13, the disk 
is much more extended and, consequently, the dwarf galaxy 
is stripped down to the stellar scales faster than, for example, 
in the reference experiment Rl. Because in all experiments 
the stars retain a certain amount of angular momentum (while 
the DM halos were non-rotating by construction) and the ori- 
entation of the angular momentum vector does not strongly 
vary with time, the stars begin to get stripped more efficiently 
compared to the DM as their angular momenta are (on aver- 
age) more strongly ali gned with th e orbital angular momen- 
tum of the dwarf (e.g. iRead et al.ll2006bl) . It is important to 
emphasize, however, that the details of how the M/L ratio 
would evolve depends, among other things, on the initial an- 
gular momentum content of the DM halos, the amount of an- 
gular momentum retained by the stars in the late stages of the 
evolution, and the degree of alignment between the angular 
momenta of stars and DM and the orbital angular momentum. 

5.4. Orbital Time or Pericentric Distance? 

In this section, we discuss the role of different orbital pa- 
rameters in influencing the transformation of disky dwarfs 
into dSphs and thus the efficiency of the tidal stirring mech- 
anism. Specifically, we focus on the relative importance of 
the orbital time, T or b, and the pericentric distance, r pe ,i, which 
determine the number and the strength of the tidal shocks, re- 
spectively. The relevant analysis is presented in Figure Q3] 
This figure shows, as a function of r per j and r o ,b, the final val- 
ues for some of the parameters that we used throughout this 
study to describe the global evolution of the dwarf galaxies 
inside the tidal field of their hosts. In particular, we consider 
the parameters r max , V max , b/a, c/a, V mt /a*, and (3. 

The filled symbols in Figure [15] present results for simula- 
tions R1-R5 which have been discussed in §[3] (see also Ta- 
ble 0). Open symbols correspond to two additional experi- 
ments, which were designed to complement the existing sim- 
ulations R1-R5. In this new set of experiments, we placed 
the dwarf galaxy model Dl used in experiments R1-R5 on the 
following orbits: (i) a short orbital time-large pericentric dis- 
tance orbit with r per i = 50 kpc, r apo = 80 kpc, and T or \, ~ 1 .7 Gyr 
and (ii) a long orbital time-small pericentric distance orbit 
with r pe ,i = 12.5 kpc, r apo = 250 kpc, and T or \, ~ 4.6 Gyr. This 
choice of parameters ensures that the pericentric distances and 
orbital times in these new simulations are not related in the 
same way as in experiments R1-R5, a fact that would be im- 
portant for the interpretation of the results. 
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FIG. 14. — Relative mass loss, \AM\/M, of the stellar (solid lines) and DM 
component (dashed lines) in simulations Rl, R2, R4, and R13. The values of 
AM/M are assigned at pericenters and are computed as \AM\/M = |(M, + [ — 
Mj)\/Mj, where M;+, and M, are the masses within r miix at the following and 
preceding apocenter, and ;' denotes the apocenter number. The relative mass 
loss of the DM initially exceeds that of the stars in all cases, reflecting the 
more efficient tidal stripping of the extended DM halos of the dwarfs. While 
the mass loss in DM dominates at all times in the reference experiment Rl, 
in simulations R2, R4 and R13 the stars begin to be stripped more effectively 
at some point during the evolution. This preferential stripping of the stars 
leads to the increase of the M/L ratio reported at intermediate and late times 
in Figuresf2]|4]and[8] 

Figure [15] confirms that the final properties of the dwarfs 
vary dramatically depending on the parameters of the orbit. 
Although the behavior of the various quantities clearly shows 
that the final outcome of the transformation depends on both 
r oi b and r per i, a closer inspection of the results in this figure 
illuminates the more crucial role of the pericentric distance 
compared to that of the orbital time in shaping dSphs via tidal 
stirring. 

Let us first consider the filled symbols in Figure [15] Dwarf 
galaxies with the same pericentric distances but with orbital 
times that differ by approximately a factor of 2 (triangle and 
pentagon) end up with a similar set of final properties. On 
the other hand, dwarfs with comparable orbital times but with 
pericentric distances that differ by a factor of 2 (pentagon, 
circle, and star) display a wide range of final properties, a fact 
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which reflects the difference in pericentric distances. 

Focusing now on the open symbols which correspond to 
the new experiments, Figure Q3] shows that the dwarf galaxy 
on the long orbital time-small pericentric distance orbit (cir- 
cle) is characterized by c/a > 0.5 and V mt /cr* < 1 and hence 
would be classified as a dSph. We note that this dwarf has 
experienced only 2 pericentric passages. On the contrary, al- 
though it has suffered a much larger number of tidal shocks, 
the dwarf on the short orbital time-large pericentric distance 
orbit (square) has c/a « 0.4 and Vrot/c* > 1 and thus did not 
undergo a transformation into a dSph. 

Overall, the results in Figure [15] lead to the following set 
of conclusions. First, the degree of the tidal transformation 
depends strongly on the combination of pericentric distance 
and orbital time. As expected, small pericentric distances and 
short orbital times, corresponding to orbits associated with a 
large number of strong tidal shocks and also characteristic of 
dwarfs being accreted by their hosts at high redshift, produce 
the most substantial tidal evolution. Furthermore, the larger 
the number of strong tidal shocks a disky dwarf suffers, the 
stronger and thus more complete is its transformation. 

Second, in order to simply transform a disky dwarf into a 
dSph (according to the criteria of § 12. 5\ without necessarily 
inducing the strongest transformation, the pericentric distance 
is a more salient parameter than the orbital time. We stress 
that this conclusion depends on the condition that T or b is short 
enough to allow the disky dwarfs to conclude at least two peri- 
centric passages 9 . Such a requirement should be generally 
satisfied in the LG environment, as a significant fraction of 
LG dwarfs have likely had enough time to complete at least 
two orbits inside their hosts. Indeed, in the context of the 
CDM cosmogony, galaxy-sized halos form early and, in par- 
ticul ar, the MW halo was probabl y already in place at z ~ 2 
(e.g.. lGovernato et al.ll2007l 12009b . In addition, structure for- 
mation in hierarchical models indicates that satellites which 
are found closer to the center of the primary system at z = 
are typically those th at fell into the host potential earlier (e.g., 
Diem and et al.l2 007). According to this, the present distances 
of most dSphs would suggest that their progenitors should 
have been accreted by their hosts at relatively early epochs. 
Assuming a MW-type host potential and the range of orbital 
parameters of LG satellites that survived until the present time 
dKlimentowski et aTll2010l) . results in a typical orbital time of 
r or b ~ 2 Gyr. Combining this value with a conservative red- 
shift of accretion of z ~ 1 suggests that dwarf galaxies in the 
LG should have been able to complete 3 or even 4 passages 
inside their hosts. This line of reasoning establishes the peri- 
centric distance, and not the orbital time, as the fundamental 
factor that controls the tidal transformation of disky dwarfs 
into dSphs in the LG and similar environments. 

5.5. Implications 

In the present study we have established that the strong 
tidal field of a MW-sized host galaxy can transform late- 
type, rotationally-supported dwarfs into stellar systems with 
the kinematic and structural properties of the classic dSphs 
in the LG and similar environments. Of course, tidal stirring 
does not constitute the only environmentally-driven mecha- 



9 As discussed in Mayer et al. 1 2001a), disky dwarfs which due to their 
very long orbital times experience just one tidal shock during their orbital 
evolution would not be transformed, even if their pericentric distances are 
fairly small (except maybe in cases of very rare, nearly radial orbits). In 
such circumstances, the orbital time will play the key role in determining the 
dynamical evolution of infalling dwarfs. 
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FIG. 15. — Final parameters of the simulated dwarfs as a function of orbital 
time (left column) and pericentric distance (right column). Filled symbols 
show results for simulations R1-R5 (see Table[2)- Open symbols correspond 
to experiments where the dwarf galaxy model Dl used in simulations Rl- 
R5 was placed on two additional orbits inside the host galaxy (see text for 
details). The pericentric distance which determines the strength of the tidal 
shocks constitutes the key factor responsible for the effective metamorphosis 
of a disky dwarf into a dSph. The combination of short orbital times and small 
pericentric distances, characteristic of dwarfs being accreted by their hosts at 
high redshift, induces the strongest and most complete transformations. 
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nism proposed for the origin of dSphs. An alternative promis- 
ing model incorporates direct gravitational interactions be- 
tween dwarf galaxies. Indeed, cosmological simulations of 
structure formation have indicated that some dwarfs may 
have en tered their Galaxy-sized hosts as members of in falling 
groups (iKravtsov et al.ll2004t ID'Onghia & Lakd f2008). It is 
thus plausible that these dwarfs could have experienced sig- 
nificant tidal pe rturbations befo r e bein g accreted by the pri- 
mary galaxies. IKravtsov et al.l (120041) did demonstrate this 
possibility in their ACDM cosmological simulations of MW- 
sized halos. 

Usin g controlled A^-body experiments, D'Onghi a et al.1 
(2009) have recently shown how the interaction between a 
pair of disky dwarfs with mass ratios ranging from 1:10 to 
1:100 can result in the rapid formation of a dSph. As the 
less massive member of the pair plunges close to the disk of 
the more massive companion, its stellar component can ex- 
perience substantial mass loss and can be tidally heated into 
a spheroidal configuration, even after only one close peri- 
centric passage. This mechanism, termed "resonant strip- 
ping", is effective when the spin angular frequency of the 
dwarf disk and the angular frequency of the orbit of the 
dwarf around the more massive system are comparable in 
magnitude, and the spin and orbital angular momenta are 
nearly aligned (that is, in a nearly prograde encounter). Such 
individual, close encounters between dwarf galaxies may 
have been more common during the assembly process of 
the host halos from acc reting groups (IKravtsov et al.l 120041 : 
ID'Onghia & Lakdl2008l) . Given that the typical velocity dis- 
persions in such groups should have been quite low (with a 
dominant group member of the size of the LMC), it is also 
possible that some of these interactions may have resulted in 
actual mergers between the dwarf galaxies. 

While all these different tidal mechanisms which rely on 
pre-processing the dwarfs before they are accreted by the pri- 
mary galaxy may play a role at some level, they must also 
be consistent with the morphology-density relation, a crucial 
constraint that the tidal stirring model can naturally satisfy. 
At present, there is no definitive observational or theoreti- 
cal evidence to rule out any of the scenarios for the origin 
of dSphs. This also includes models that advocate the ab ini- 
tio for mation of dSphs from cosmological initia l conditions 
(e.g- lRicotfi & Gnedinll2"00l ISawala et al.ll2010l) . In all like- 
lihood, all mechanisms do operate simultaneously to differing 
degrees in producing the population of dSphs in the LG and 
similar environments. 

Our investigation of the role of different pericentric dis- 
tances and orbital times on the efficiency of tidal stirring has 
interesting implications for the expected properties of dSphs. 
Our results show that in order to transform a disky dwarf with 
an orbital apocenter in excess of 200 kpc into a dSph, an orbit 
with an eccentricity that is substantially larger than the me- 
dian expected eccentricity is required (see Figures l2l and n~5T>. 
Given that the current distances of dwarf galaxies in the LG 
are suggestive of their orbital apocenters, the distant dSph Leo 
I (currently at a di stance of ~ 250 kpc from the MW; e.g., 
ICaputo et al.lll999r) should be on a very wide, nearly radial 
orbit, assuming that it is the product of tidal stirring. Fur- 
thermore, even on such eccentric orbits the remnants possess 
some residual rotation at the level of Vrot/c* ~ 0.5 (see open 
circle in Figure [L5"ll. We thus predict that Leo I should have a 
higher V ro t/o'» compared to the dSphs located closer to the pri- 
maries. Interestingly, albeit in conclusive, there are signatures 
of intrinsic rotation in Leo I (Soh n et al.ll2007t iMateo et"aT] 



2008). We note that such rotation has b een shown to exp lain 
the observed kinematics of this dwarf (Lokas et al. 2008). 

Our prediction that in the context of tidal stirring wider or- 
bits should be associated with less evolved dSphs, and thus 
with higher values of Vrot/o"*, should apply, even more so, to 
the oddly isolated galaxies Cetus and Tucana, which are lo- 
cated in the outskirts of the LG at more than 500 kpc from 
M31 and the MW, respectively. Although measurements of 
kinematics are extremely challenging due to the very large 
distances involved, recent spectroscopic studie s have pre- 
sented tentative ev idence for rotation in Cetus (iLewis et al.l 
120071) and Tucana dFraternali et al.l l2009) with Vrot/o - * values 
in the range of 0.5 < Vi-ot/c* < 1. If confirmed by future ob- 
servations, detection of rotation at such levels would be con- 
sistent with the predictions of the present study. 

Furthermore, at distances in the range ~ 300 - 600 kpc 
from the MW and M31, we find the dwarf galaxies Phoenix 
and LGS3, w hich are classified as "transition-type" dwarfs 
(Mateo 1998). While these objects exhibit structural proper- 
ties consistent with those of d Sphs, they are associated with 
some residual gas (see, e.g., iGrcevich & Putmanl 120091 and 
references therein). On the other hand, these transition-type 
dwarfs exhibit negligible or no ongoing star formation (no de- 
tection of Hjj regions has been reported), a fact that differ- 
entiates them from standard gas-rich, star-forming dlrrs (e.g., 
iGrebel 1999; Skillman 2005). Being at intermediate distances 
between dSphs and dlrrs and combining properties of both 
classes of objects, these transition-type dwarfs can provide 
valuable tests for the predictions of tidal stirring regarding the 
different stages of the transformation process. 

One possibility for their origin is that their progenitors were 
disky satellites that were accreted only recently by their hosts 
on wide, fairly radial orbits, and are still in the process of 
being transformed into dSphs by tidal stirring. It is then plau- 
sible that these transition-type dwarfs have only concluded 
one pericentric passage (or they are on the way to their first 
pericentric approach), having lost most of their gas by ram 
pressure stripping. This is a reasonable assumption given that 
gas removal by the combination of tides and ram pressure pro- 
ceeds faster tha n the morphologica l transformation of the stel- 
lar component dMayer et alJ l2006). If the scenario proposed 
here is correct and their transformation has only been partial, 
transition-type dwarfs should possess significant rotation in 
their stellar component (as in the case of the remote dSphs 
Leo I, Tucana, and Cetus discussed above). Verifying obser- 
vationally this basic theoretical prediction would be impor- 
tant - unfortunately, the stellar kinematics of these systems is 
poorly known. Moreover, provided that they are at the inter- 
mediate stages of their transformation, our results suggest that 
they should also show signs of bar-like distortions. Definitive 
conclusions on all of these issues would require a synergy be- 
tween detailed photometric and kinematic measurements as 
well as exhaustive comparisons with theoretical models. To 
this end, extending the search for evidence of tidal stirring 
in systems with more recent assembly histories than the LG 
may reveal a population of satellites in the process of being 
transformed, and thus offer unique opportunities to constrain 
the tidal stirring model. Interestingly, recent studies using the 
SDSS database have uncovered a population of dwarf galax- 
ies in the Virgo cluster that exhibit disk-like features and bars, 
appearing t o be in a tr ansitional stage between a disk and a 
spheroid dLisker et alJl2007h . The existence of this class of 
objects is, in broad terms, consistent with the predictions of 
tidal stirring for the different stages of the transformation. 
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Several alternative scenarios to tidal stirring exist for ex- 
plaining the nature of the spheroidal morphology of the dis- 
tant LG dwarfs. For example, as a result of three-body 
interactions, satellites can acquire extremely energetic or- 
bits with apocenters beyond the v irial radius of the primary 
and be ejected to large distances (Sale s et aT1l2007l) . In this 
model, ejected subhalos are typically the less massive mem- 
bers of a pair of satellites that is tidally disrupted during the 
first approach onto the host. Such ejections can also oc- 
cur during the tidal disruption of a bound system of mul- 
tiple subhalo s that is accreted a s a single unit by the pri- 
mary galaxy (Ludlow et al. 2009). In both of these scenarios, 
strong tidal interactions with the more massive companions 
can took place before the satellites are accreted by their hosts. 
Such encounters may eventually tr ansform the dwarfs i nto 
dSphs, either via resonant stripping (D'Onghia et al. 2009) or 
other gravitational pr ocesses such as mergers. Interestingly, 
iKravtsov et alj d2004l) reported in their cosmological simula- 
tions the existence of a few satellite systems at distances of 
~ 1000 kpc from their MW-sized hosts. Based on their fi- 
nal Viot/er» values, these objects would be classified as dSphs 
and could plausibly represent the counterparts of Cetus and 
Tucana in the context of the LG. The tidal heating of these 
systems did occur in small groups that were accreted by the 
primary halo at the present epoch. As the tidal forces unbind 
these accreting groups, energy redistribution can increase the 
orbital energy and apocentric distance for some of the satel- 
lites, providing an explanation for the presence of isolated 
dSphs in the periphery of the host galaxy . The extreme radial 
velocity of Leo I (e.g., iMateo etaTll2008l) and the rela tively 
high recession velocity of Tucana dFraternali et al .1120091) may 
already suggest that these systems have been propelled into 
their highly energetic orbits th roug h the type of intera ctions 
suggested by lSales et aD (120071) and lLudlow et alJ (120091) . 

Binary mergers between individual satellite galaxies tak- 
ing place outside of infalling groups may offer another alter- 
native explanation for the puzzling presence of the isolated 
dSphs in the LG. Indeed, it has already been demonstrated 
that interactions and direct merge rs of subhalos can l ead to 
their very strong evoluti on (e.g., iKnebe et al. 2004, 120061: 
Klime ntowski et aO2010l) . For example. iKlimentowski et alJ 
(2010) found in their constrained simulation of the LG that 
~ 10% of all surviving subhalos in the MW and M3 1 have un- 
dergone a major encounter with another subhalo in the past. 
Most of these events occurred at very early times, between 
Z ~ 1 and z ~ 3, while the interacting subhalos have not yet 
become satellites and are still outside their hosts. Whether 
binary mergers between dwarfs can explain the existence of 
dSphs orbiting far from the primary galaxies, such as Tu- 
cana and Cetus, is currently under investigation via a com- 
bination of cosmological simulations and controlled numer- 
ical experiments (Kazantzidis et al., 2010a, in preparation). 
Ongoing investigations aimed at studying in detail the stel- 
lar components of these isolated dSphs (e.g., [Bernard et al] 
l2009tlMonelli et alJ2010al). includ ing their star formation his- 
tories (e.g.. iMonelli et aLll2010bl) . may soon provide useful 
constraints on the various competing models that have been 
proposed so far for their origin. 

In this paper, we have investigated the tidal evolution of 
dwarf galaxies comprising exponential stellar disks. Recent 
cosmological simulations do support the idea th at isolated 
dwarf s are rotating disk s dGovernato et al.l 120101 : see, how- 
ever, [Sawajaetay|20T(J). The morphological and dynamical 
transformation of such disky dwarf galaxies into dSphs under 



the action of tidal forces is a rather rich process which in- 
volves several stages and is characterized by different events. 
Consequently, numerical simulations where a sp heroidal stel- 
lar system is postulated from the beginning (e.g., Mun oz et al.l 
2008; Pen arrubia et al.ll2008l) may be inadequate to describe 
how dSphs have evolved to the present time subject to the 
tidal field of their hosts. Indeed, if the picture presented here 
is correct, dSph galaxies should have formed relatively late in 
most cases and only after concluding a number of pericentric 
passages inside the primary galaxies. 

More specifically, Penarrubi a et al.l d2008) performed a se- 
ries of A^-body simulations to study the dynamical evolu- 
tion of dSphs in a host potential assuming a spherical King 
model for the stellar component embedded within an NFW 
halo. These authors reported that tidal effects lead to an in- 
crease in the M/L ratio in most cases, in contrast to the find- 
ings of the present study. This discrepancy is possibly due 
to the fact that the stellar distribution in the Penarrub ia et al.l 
(2008) dwarfs followed the King profile which is character- 
ized by a density core. As a result, the stars were loosely 
bound within the potential of the dwarf galaxy and the stel- 
lar component was much more heavily stripped compared to 
what is typically found here. In our experiments, the forma- 
tio n of tidally-induced bars , which is obviously missing from 
the Penarrubia et al. (2008) models, is crucial. This is because 
bar formation enhances the resilience of the dwarf galaxies to 
mass loss and tidal stripping by increasing the stellar density 
and, correspondingly, the depth of the potential well. 

Given the existence of a number of qualitatively different 
scenarios for the formation and evolution of dSphs, it is crit- 
ical to be able to test these scenarios and discriminate among 
them. The growing kinematic data sets for dSph galaxies 
can already facilitate detailed c omparisons with theoretical 
models (e.g., Walker et al. 2009). For example, the level of 
residual stellar rotation in the remnants can be used to con- 
strain the competing models. Indeed, according to the re- 
sults presented i n Figure [TT1 and de tailed analysis of similar 
simulations (see lLokas et alj|2010al) . if the present-day dSph 
galaxies originated from disky dwarfs they should, at least 
in some cases, show signatures of intrinsic rotation. Inter- 
es tingly, in addition to the cases of th e isolated dSphs Leo 
I (So hn et all 12007b IMateo et al.l 12 008). Cetus (iLewis et al l 
120071) . and Tucana (IFraternali et al.. 2009) mentioned already, 
detection of rotation at different le vels of significance has 
also been claimed for Ursa Minor dHargreaves et al.l 1 1 994t 
I Armandroff et all 1 1 9951) and Sculptor (Battaglia et al.M2Q08l) . 
It is important to stress, however, that intrinsic rotation may 
be difficult to distinguish from the velocity gradients induced 
either by the presence of tidal tails or, for nearby systems, 
caused by t ransverse motions ( the so-called "perspective ro- 
tation"; see IWalker et al.ll2008l) . Due to such complications, 
Leo I seems to be the most promising c andidate for detection 
of intrinsic rotation (see discussion in lLokas et al.ll2008l) . In 
addition, as discussed above, our simulations predict a pos- 
itive correlation between the magnitude of stellar rotation in 
dSphs and their distances from the center of the host galaxy 
(or equivalently the time of infall of the progenitor dwarf onto 
the primary). Such a correlation would be difficult to estab- 
lish, for example, in models that propose the ab initio for- 
matio n of isolated dSphs from cosmological initia l conditions 
(e.g.. lRicotti & Gnedinll2005l: ISawala et : al.ll2Q10t) . Indeed, in 
this case the amount of rotation in the stellar components of 
dSphs should depend on the specifics of the formation process 
and the intrinsic properties of these systems (e.g., the initial 
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spin parameter of the DM halo of the dwarf) and thus be un- 
correlated with the distance from the center of the host galaxy. 

In the context of specific observational signatures of the 
tidal stirring model, bars play a prominent role. According 
to our results, bar-like structures should be relatively com- 
mon in dSphs as the bar phase is one of the longest stages in 
the transformation process. In fact, some of the less evolved 
dSphs in the LG may still be in the bar stage and show signs 
of bar-like distortions. One irrefutable example of a LG dwarf 
that contains a bar is the LMC. Orbital evolution models us- 
ing three-dimensional v elocities constrain ed by recent proper 
motion measurements dKallivayalil et alJl2006l) suggest that 
the LMC may currently be on the first passage having just 
cross ed the pericenter of its orbit around the MW (Bes la et alj 
2007). Confirming that the LMC is on a such extended or- 
bit with a pericentric distance of r per i ~ 50 kpc would be in 
line with our findings, as we expect its stellar component to 
be less evolved and thus in some early, transitory stage of the 
transformation between a disk and a spheroid. 

Except for the LMC, however, the overall number of de- 
tections of bar-like structures among the dSphs in the LG 
is relatively low 10 . Indeed, markedly elongated isophotes 
that could be attributed to a residual bar-like component 
are observed in only a few cases, including Ursa Minor 
dlrwin & Ha tzidimitriou 1995) and the newly discovered Her- 
cules dSph (IColeman et al.ll20071) . Recently, the strongly non- 
spherical shape of the core of the disrupting nearb y Sagittarius 
dSph h as also been ascribed to a pre-existing bar (Lokas et al. 
2010b). Interestingly, the recently discovered ultra-faint MW 
satellites als o exhibit substantial degrees of flattening (e.g., 
MaitmetaL||2008). Of course, the high elongation of the 
stellar component in some of the dSphs can also be caused 
by other effect s, including the triaxia lity of the surround- 
ing DM halos (Kaz antzidis et al.l 120 1 Oh or tidal deformation 
in the gravit a tional field of the MW. Regarding the latter, 
Martin & Jin (2010) have recently proposed that Hercules 
dSph is in fact a stellar stream in formation, thus suggest- 
ing tidal disruption as the most valid scenario for the extreme 
shape of this system. Deeper observations to track evidence 
of such tidal interactions would be able to settle these issues. 

The findings of the present study also indicate that 
the formation of tidally-induced bars is strongly linked 
with the transformation of rotationa lly-supported dwarf s to 
dSphs (for similar conclu sions, see iMayer et al.l 1200 lal and 
iKlimentowski et ai1l2009al) . Indeed, as shown in Table |3] in 
only one of the 14 simulations that produced dSphs in the end, 
specifically experiment R16, a bar did not form at some point 
during the evolution of the progenitor disky dwarf inside the 
host galaxy. The following sequence of events is typically 
observed in our simulations. First, the strong tidal forces at 
the initial pericentric approach trigger a bar instability in the 
disk of the dwarf. The tidally-induced bar transports angu- 
lar momentum outwards to the outer regions of the disk and 
to the DM halo. As tidal stripping removes the outer parts of 
the dwarf, the entire angular momentum content gradually de- 
creases and the ability of the dwarf galaxy to be supported by 

10 We note that there are a number of reasons as to why bars in dwarf galax- 
ies may escape detection in the LG. First, bars are oriented randomly with re- 
spect to the observer, so some bar-like structures may appear as only slightly 
flattened. Second, the smoothing procedures usually applied when measur- 
ing the surface density distributions of stars in dwarf galaxies can decrease 
the detectability of bars in these systems (for a thoroug h discussion pertain- 
ing to the difficulties in identifying bars in dSphs, see Klimentowski et al. 
2009a). 



rotation progressively diminishes. Second, subsequent tidal 
shocks destroy the centrophilic stellar orbits which support 
the bar and increase the stellar velocity dispersion. As a re- 
sult, the bar continuously loses its elongation and is tidally 
heated into a more isotropic, diffuse spheroid. The ultimate 
outcome of these two physical processes is the formation of 
pressure-supported stellar systems with values of V mt /a* < 1 
that are appropriate for dSphs. 

We stress that the above discussion describes only one 
channel for the transformation of a disky dwarf into a dSph via 
tidal stirring. An alternative mechanism involves the buckling 
of the bar d ue to the amplificat ion of vertical m = 2 (bend- 
ing) modes (IMayer et al.ll2001allbl) . While the loss of angu- 
lar momentum occurs in exactly the same way as previously 
described, the buckling of the bar now becomes the main 
process that leads to the significant vertical heating and the 
increase of the veloc ity dispersion of the system (see, e.g., 
Debattista et al. 2006 for a detailed description of these pro- 
cesses). In this picture, the spheroidal shape of the stellar 
distribution constitutes the end result of the instability and is 
not driven by tidal heating as in the simulations of the present 
study. Because the buckling instability requires a fairly strong 
bar to develop, this transformation mechanism is relevant to 
relatively massive systems with high surface densities. In- 
deed, none of the bars in our low surface density dwarfs 
showed any signs of buckling. 

To summarize, it seems that there are at least two channels 
for the formation of the spheroidal component in the tidal stir- 
ring scenario: prolonged impulsive tidal heating and the bar 
buckling instability. The former is favored in low mass, low 
surface density disky dwarfs for which heating is particularly 
efficient. This mechanism likely applies to the progenitors of 
the faintest classic dSphs such as Draco, Sculptor or Leo I. On 
the other hand, bar buckling requires strongly self-gravitating 
systems with higher surface densities. This channel plausi- 
bly operates on the progenitors of the brightest dSphs such 
as Fornax and Sagittarius, and even more so on those of the 
bright dwarf elliptical satellites of M31, such as NGC185 and 
NGC167. This mechanism was also found to be applicable 
to the transformation of re latively massive and brig ht spiral 
galaxies in galaxy clusters (Mastropietro et al. 2005|) 

Our simulations also suggest that the amount of mass loss 
that the dwarf galaxies experience can be considerable. In- 
deed, in the reference experiment Rl, the dwarf lost ~ 90% 
of its initial mass within r max and still survived as a bound 
entity (the amount of mass loss for the most heavily stripped 
dwarfs in simulations R2 and R4 reached ~ 99%). Typically, 
the maximum circular velocities, V max , decreased by a factor 
of ~ 2 during the orbital evolution in the cases where dSphs 
were produced (this factor increased to ~ 3 in experiments 
R2 and R4). Such substantial evolution in V max occurred de- 
spite the presence of the baryons which tend to moderate the 
effect of tidal shocks by making the potential well deeper, es- 
pecially after bar formation. These findin gs are in agreement 
with those of Klime ntowski et all (l2009al) as well as with re - 
sults of other studies. For example. iRravtsov et al.1 (120041) . 
using fully cosmological DM-only simulations, reported that 
the V max of their most evolved subhalos decreased by a factor 
of 1 .5 - 2 on average during 1 G yr of tidal evolution . More 
recentlv. lDiemand et al]d2007l) and lMadau et al.l(l2008l) found 
that satellites in the Via Lactea simulation lost between ~ 30 
and ~ 99% of their pre-infall mass, and that V max typically 
decreased by a factor of ~ 2 - 3. 

The previous discussion suggests that for subhalos which 
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have been accreted by their hosts at early cosmic epochs and 
have completed several orbits with fairly small pericenters, 
Vmax i s expected to have evolved significantly. Therefore, 
the DM halos of present-day dSphs may have had consid- 
erably larger masses and circular velo cities when they en- 
tere d the halo of the M W (see, e.g., iKravtsov et alJ 120041 
and iMadau et al.l 120081) . This has important implications 
for galaxy formation models. Indeed, the substantial obser- 
vational work during the past few years suggests that the 
present-day classic dSphs of the MW and M3 1 have central 
stellar vel ocity dispersions in the range c* ~ 7- 13 kms" 1 
(see, e.g., IWalker et al.1 12009 and references therein). As- 
suming that V max ~ v3(7*, in accordance with the results 
of § 15.11 implies current values of V max in the range V max ~ 
12-22 kms" 1 and initial V max values, namely before infall 
and tidal mass loss, in the range V max ^24-44 kms" 1 . These 
numbers are quite important. Indeed, with such relatively 
high values of initial halo V max , photoevaporation of the gas by 
the cosmic ultraviolet background after reionization as well as 
supernovae feedback could not be the factors that shaped the 
baryonic content and nature of these dSph s. In fact, accord- 
ing to the radiative transfer simulations of ISusa & Ume mural 
(2004), photoevaporation would be effective and remove most 
of the gas only for V max < 20 kms" 1 . This is consistent with 
the lack of a clear signature of th e reionization epoch in the 
star formation h istories of dSphs (Grebel & Gallagherj |2004t 
lOrban et al.ll2008l) . 

In this paper, we have also demonstrated how the orbital pa- 
rameters and initial structures of the progenitorlate-type disky 
dwarfs determine the final properties of the dSphs. There- 
fore, the fact that Fornax and Draco have roughly similar 
masse s at present, as inferred from the V max of their halos 
(e.g.. iKazantzidis et al.l 1200 4b). but differ by about an order 
of magnitude in luminosity and M/L ratio, can be explained 
in two ways. One possibility is that the progenitors of these 
two dwarfs began with very different relative amounts of DM 
and baryons, for reasons related to their formation history and 
not to the environment. Alternatively, Fornax and Draco orig- 
inated from systems that had comparable DM and baryonic 
masses, but experienced dissimilar tidal evolutions because 
they entered the primary galaxy at different epochs and/or on 
different orbits. 

Regardin g the latter, accordi ng to the hydrodynamical sim- 
ulations of Ma yer etail ((2007), gas-rich disky satellites that 
were accreted by their hosts when the intensity of the cosmic 
ultraviolet background was much higher than today (z ~ 2), 
can be completely stripped of their gas by ram pressure in 
one or two pericentric passsages. In this case, the final sys- 
tems would correspond to dSp hs with truncated st ar formation 
histories such as Draco ("e.g.. lOrban et a l. 2008). While the 
baryonic content of the progenitor dwarfs decreased signifi- 
cantly as a result of gas stripping, the initial DM mass in the 
central regions around the surviving baryonic core is largely 
preserved. This is because DM is affected only by tides and 
not by ram pressure. As a result, exceptional M/L values, 
of the order of 100, similar to those inferred for Draco and 
Ursa Minor, can arise. On the other hand, disky dwarfs that 
were accreted by their hosts at z < 1, when the intensity of the 
cosmic ultraviolet radiation dropped by more than an order of 
magnitude compared to z ~ 2, were able to retain some gas be- 
cause tides and ram pressure could not strip it completely. Un- 
der these conditions, the infalling dwarf galaxies can undergo 
tidally-triggered bursts of star formation associated with bar- 



driven gas inflows at pericentric approaches. Such a model 
would be applicable to dSphs with extended star formation 
histories such as Fornax, Carina, and Leo I. This mechanism 
can produce dSphs that are brighter for a given halo mass (or 
central stellar velocity dispersion) compared to the ones that 
were accreted earlier. While this generic scenario seem to ex- 
plain the differences in the properties of Fornax and Draco 
despite their similar masses, it would be important to revisit 
it in future work with models capable of capturing the multi- 
phase structure of the ISM in dwarf galaxies. 

Lastly, it has long been debated whether the inferred high 
M/L ratios of dSphs indeed signify exceptional DM content 
or are simply a reflection of strong tidal effects and of the fact 
that these systems are in reality on the verge of disruption. 
The findings of the present study have intriguing implica- 
tions for earlier attempts to model dSph galaxies as unbound 
system s without DM (e.g., Kroupa] [l997t iKlessen & Kroupal 
1998). Indeed, the tidal stirring model demonstrates how 
mass loss and the formation of tidal tails can be consis- 
tent with the presence of a bound stellar component em- 
bedded in a relatively massive CDM halo, even after sev- 
eral Gyr of tidal evolution inside the host. Our results indi- 
cate that substantial mass loss and the existence of a gravita- 
tionally bound dSph galaxy with a relatively high M/L ratio 
are not mutually exclusive, confirming earlier claims based 
on a s maller set of lower resolution simulations (May er et al.1 
2002). The findings of the present study also suggest that the 
claimed detection of extra-tidal stars in a number o f dSphs, 
including Ursa Minor (Martmez-Delgado et al. 2001), Fornax 
dColemanetal] 120051) . Carina (iMufioz et al.l l2006h . and Leo 
I (ISohn et al.1 120071) is consistent with satellite accretion in 
CDM models. We note in passing that the number of such 
detections may be low due to the intrinsic difficulties associ- 
ated with separati ng the tidal tails from the bo und core of the 
dwarfs (see, e.g., iKlimentows ki et al.l l2009b). Furthermore, 
in the context of tidal stirring, dSphs embedded in CDM ha- 
los exhibit stellar distributions that are adequately fit by ex- 
ponential or Ki ng profiles in agreement with observed dSphs 
(see review by Maverl l2010l) . Stellar profiles of this type are 
difficult to accommodate within models that represent dSphs 
as systems devoid of DM. Indeed, such models predict nearly 
flat profile s as expected for objects close to complete disrup- 
tion (e 



mles as expected 1 
.g., Kroupa 1997). 



5.6. Caveats and Future Directions 

Certainly the approach presented in this paper is not with- 
out caveats. A first limitation is related to the fact that we 
have adopted a single primary galaxy with the present-day 
structural properties of the MW. In general, at the time when 
the dwarfs are accreted by the primary galaxy at high red- 
shift, the DM and baryonic masses of the host would be dif- 
ferent compared to those of the present time Our methodology 
thus neglects the cosmological evolution of the host galaxy 
structure via mergers and smooth accretion during the inter- 
actions with the dwarfs. However, this simplification may 
be justified to a certain degree by recent ACDM galaxy for- 
mation simulations, which have shown that MW-sized disk 
galaxies assemble most of the mass in their inner regions be - 
tween - 8 - 10 Gyr ago (e.g., iGovernato etail 120071 120091) . 
This is indeed the timescale that we follow in our simulations. 
Nonetheless, a more complete investigation would have to in- 
clude the ongoing formation of the host galaxy. In addition, 
we have assumed a host DM halo that is spherical (except in 
the very inner regions that are dominated by the potential of 
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the di sk), instead of tri axial as postulated by CDM models 
(e.g.. iFrenk et al.lll988l) . Halo triaxiality and the complex- 
ity of halo formation in a realistic cosmological context, with 
continuous mergers, accretion, and rapidly changing potential 
wells can have an impact on the orb ital evolution of infalling 
satellites (e.g.. iRravtsov et aI1l2004l) . with consequences for 
the efficiency of their transformation. It will be important to 
explore these issues in future investigations of tidal stirring. 

A second shortcoming of our work is related to the fact that 
we have focused on experiments where the alignments be- 
tween the internal angular momenta of the dwarfs, those of the 
primary disks and the orbital angular momenta were all pro- 
grade (45° in most cases; see Table O. Prograde alignments 
between the orbital angular momenta of the dwarfs and the 
spins of the primary disks, or between the angular momenta 
of the two disks, are expected to be important mainly for or- 
bits with pericentic distances smaller than a characteristic ra- 
dius containing a non-negligible fraction of the mass of the 
primary disks. In this case, the tidal effects of the host disks 
would be enhanced and the orbits of the dwarfs would quickly 
decay due to the additional dynami cal friction provided by th e 
disks of the primary galaxies (e.g., Ouinn & Goodman 1986). 
As our typical pericentric distances are r per i > 15 kpc, we ex- 
pect the efficiency of the transformation reported here to be 
weakly affected by such effects. 

On the other hand, the orientation between the orbital an- 
gular momenta of the dwarfs and the internal angular mo- 
menta of their disks is particularly relevant for our experi- 
ments. For retrograde alignments between these two angu- 
lar momenta, tida l stripping is considerably reduced (e.g., 
iRead et aill2006bl) and both bar formation and the efficiency 
of transfor mation via the tidal stirring mechanism are also 
suppressed (Mayer et al. 2001a). These facts highlight the 
importance of coupling between orbital and internal mo- 
tions. Determining the statistics of alignments between in- 
ternal and orbital angular momenta would require a series 
of fully cosmological studies, focusing on the accretion of 
satellites within halos of disk galaxies, and is therefore be- 
yond the scope of this paper. We stress, however, that the 
orbital parameters adopted in the present study correspond to 
only moderate alignments, and therefore our results should 
not be biased by any strong coupling of angular momenta. 
Unless there is some yet unknown cosmological bias against 
the mildly prograde alignments that we have adopted here, 
our findings should be able to illuminate at least some of the 
details of the typical transformation experienced by infalling 
disky dwarfs in the LG and similar environments. 

Lastly, the most evident limitation of this study is that we 
have addressed the efficiency of the tidal stirring mechanism 
only in the collisionless regime. Our results should therefore 
be viewed as preliminary. A more complete treatment includ- 
ing hydrodynamics is required to illuminate one of the most 
distinct pr operties of dSphs, name ly their low gas content 
(see, e.g., iGrcev ich & Putman 2009 and references therein) 
and fully refine the conclusions presented here. The present- 
day structure of dwarf galaxies originates from a complex 
interplay of effects and a full explanation requires detailed 
knowledge of their star formation histories and chemical evo- 
lution, amongst others. Adding star formation as a further in- 
gredient will offer the possibility to determine the magnitude 
of starbursts induced in the dwarfs at pericentric passages, 
while gaining a deeper understandi ng of the wide diversity in 
their star formation histories (e.g., Grebel 2000; Orban et al. 
2008). Furthermore, specific predictions for the metallic- 



ity of the remnants formed by tidal stirring, let alone com- 
parisons with t he luminosity-metal licity relation for nearby 
dwarf galaxies dTolstov et al.ll2009l) . are not currently avail- 
able. We plan to extend the present investigation in these di- 
re ctions in due course. 

May er etafl J20 07) have also demonstrated how the effi- 
ciency of tidal stirring is affected by the presence of a dis- 
sipative component, and how it varies depending on the bal- 
ance between heating and cooling in the gas. This is espe- 
cially relevant for disk-like progenitors that are able to over- 
come ram pressure and retain their gas for a longer period of 
time. Satellites infalling at z < 1 when the intensity of the 
cosmic ultraviolet background is weaker and the gas can set- 
tle in a colder and denser configuration within the potential of 
the dwarf would fall in this category. In this case, the tidal 
heating of the bar into a spheroid becomes less efficient. This 
is because the bar instability causes the gas to flow towards 
the central region of the dwarf, increasing its central density 
and causing a more adiabatic response of the system to the 
external tidal perturbation. 

Despite the aforementioned limitations, several facts do 
suggest that our results for the efficiency of tidal stirring 
should be regarded as conservative. First, as we discussed 
above, some of the progenitors of present-day dSphs might 
have suffered significant ti dal perturbations befo re being ac- 
creted by their hosts (e.g.. IRravtsov et aLll2004h . Therefore, 
it is plausible that these systems might have entered their pri- 
mary galaxies already partially transformed. Such a condition 
would facilitate their complete transformation inside the tidal 
field of the primary galaxy. Second, if the progenitor disky 
dwarfs had either very low halo concentration parameters or 
core-like density profiles, as it is indeed suggested by both 
the modeling of rotation cur ves of present-day LSB and din- 
galaxies (e.g. Jde Blokll2010l) and recent cosmological si mula- 
tions of dwarf galaxy formation dGovernato et al.ll2.010t) , our 
conclusions regarding the effectiveness of the transformation 
would be reinforced. Indeed, in such circumstances the re- 
sponse of the dwarfs to the tidal shocks would be much more 
impulsive compared to the cases of steep NFW-like profiles 
like the ones that we used in our experime nts. This would 
give ri se to augmented tidal mass loss (e.g., Kazantzidi s et al.1 
2004a) and to a more effective transformation into a dSph. 
Lastly, none of the disky dwarfs in our experiments passed 
through the disk of the host galaxy. This is important as it has 
been recently shown that disk shocking, namely tidal shocks 
induced by passages through the disk, can affect significantly 
the evolution of satellites having masses < 1O 9 M and peri- 
centic distances < 30 kpc, and even ca use the disruption of a 
fraction of them (lD'Onghia et al.ll2010l) . 

As a final remark, we reiterate that accumulating observa- 
tional and theoretical evidence suggests that dwarf galaxies 
are not formed as t hin disks, but rather a r e born as thick, 
puffy systems (e.g., iDalcanton et al.l l2004t Kaufman n et al.l 
120071) . The effect of thermal support, as opposed to rota- 
tion, ought to be thoroughly investigated in forthcoming stud- 
ies of tidal stirring. Indeed, a thicker, more diffuse stel- 
lar c omponent sug gests a stronger effect of direct tidal heat- 
ing (ISpitzerl Tl958h . but the bar instability, an essential ele- 
ment of tidal stirring, is associated with thin stellar config- 
urations. Moreover, gas stripping by ram pressure should 
be enhanced in an initially thicker, more diffuse stellar sys- 
tem owing to the reduced gravitational restoring force of the 
gas. The interplay between all these aspects of the model- 
ing will be assessed with future work where realistic gas-rich 
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dlrrs will be evolved inside MW-sized hosts with recipes of 
radiative cooling, star formation, and supernovae feedback 
(Kazantzidis et al., 2010b, in preparation). These dwarfs have 
been formed sel f-consistently in cosmological hydrodynami- 
cal simulations (iGovernato et al.ll2010l) and are characterized 
by a turbulent and multi-phase realistic ISM. 

6. SUMMARY 

Using a suite of collisionless A^-body simulations we have 
investigated the efficiency of the tidal stirring mechanism for 
the origin of dSphs. Specifically, we have examined the de- 
gree to which the sizes, masses, shapes, and kinematics of 
late-type, rotationally-supported dwarfs are affected by the 
gravitational field of MW-sized host galaxies for a range 
of dwarf orbital and structural parameters. Unlike previous 
work on the subject, we have employed equilibrium numer- 
ical models of dwarf galaxies constructed from composite 
DFs and consisting of exponential stellar disks embedded in 
cosmologically-motivated DM halos. The self-consistency 
of the adopted models is crucial for confirming the complex 
transformation process of a disky dwarf into a dSph. This 
aspect of the modeling constitutes the major improvement 
that we introduce in the present study. Furthermore, we have 
extended earlier contributions on the subject by conducting 
a simulation campaign which is carefully designed to allow 
an investigation of a much larger parameter space than be- 
fore. Lastly, the fairly high numerical resolution of our ex- 
periments combined with the growing observational data sets 
for dSph galaxies provide unique opportunities for systematic 
and quantitative comparisons between the theoretical models 
and the data, and we undertake such a task in a companion 
paper (Lokas et al. 2010, in preparation). 

Our main results and conclusions can be summarized as fol- 
lows. 

1. Tidal interactions between rotationally-supported 
dwarf galaxies and MW-sized hosts can lead to the 
formation of pressure-supported, spheroidal stellar 
systems with kinematic and structural properties akin 
to those of the classic dSphs in the LG and similar 
environments. Our exploration of a wide variety of 
initial conditions for the progenitor disky dwarfs sug- 
gests that such tidal transformations are fairly efficient 
and should thus be common occurrences within the 
currently favored cosmological paradigm. Due to 
the fact that satellite accretion is a generic feature 
of hierarchical models of structure formation, the 
transformation process described in this study should 
be applicable to at least some of the dSph galaxies in 
the universe. 

2. The transformation mechanism is complex and involves 
a combination of tidally-induced bar instabilities in 
stellar disks and impulsive tidal heating of the stellar 
distribution. Given the self-consistency of our dwarf 
galaxy models, we can safely conclude that the forma- 
tion of dSphs can be entirely attributed to tidal pertur- 
bations, rather than being a consequence of the initial 
conditions. In the context of the tidal stirring model, bar 
formation is intimately linked to the formation of dSphs 
(see Table©. Loss of angular momentum caused by the 
bar instability and simultaneous increase of the stellar 
velocity dispersion due to tidal heating lead to low val- 
ues of Vrot/c, in the simulated remnants comparable 



to those of observed dSphs (Vrot/c* ^1)- Heating via 
tidal shocks at pericentric passages decreases continu- 
ously the elongation of the bar and causes the initially 
disky stellar distributions to transform into spheroids 
with projected axis ratios of c/a > 0.5. 

3. Bar formation constitutes a sufficient but not necessary 
condition for the formation of dSphs via the tidal stir- 
ring of rotationally-supported dwarfs. In cases where 
bar instabilities are not triggered by the tidal interac- 
tions with the host galaxies, spheroidal stellar systems 
with negligible amounts of rotation can still be pro- 
duced solely via the action of impulsive tidal heating 
(R16; see § |4T2l 

4. The effectiveness of the transformation into a dSph 
depends crucially on the orbital parameters of the 
progenitor disky dwarfs (§ (3}. For a fixed eccen- 
tricity, r a p /r pel -i, tighter orbits which are character- 
ized by shorter orbital times, r or b, and smaller apoc- 
entric, r apo , and pericentric distances, r per i, lead to 
more rapid and complete transformations. For a 
fixed apocentric distance, orbits with higher eccentric- 
ities also induce stronger transformations. Under the 
right combination of orbital parameters, tidal stirring 
can yield spherically-symmetric (b/a w c/a > 0.95) 
and isotropic (J3 m 0) stellar systems with negligible 
amounts of rotation (Vrot/c* ^ 0.1) (R2 and R4; see 
Table©. 

5. The degree of the tidal transformation depends on both 
the number and the strength of the tidal shocks, which 
are determined by the orbital time, T or b, and the pericen- 
tric distance, r per ;, respectively. Small pericentric dis- 
tances and short orbital times, corresponding to orbits 
associated with a large number of strong tidal shocks, 
produce the strongest and most complete transforma- 
tions. However, in order to simply transform a disky 
dwarf galaxy into a dSph without necessarily inducing 
the strongest transformation, r pei i is the more salient or- 
bital parameter. The last conclusion holds provided that 
r or b is short enough to allow the dwarfs to complete at 
least two pericentric passages inside their hosts. Such a 
condition is likely satisfied for a significant fraction of 
dwarf galaxies in the LG (see § 15.4b . 

6. The efficiency of the transformation via tidal stirring 
is notably affected by the structure of the progenitor 
disky dwarfs. Specifically, it is enhanced considerably 
for those with less massive and more extended disks, as 
well as for dwarfs embedded in halos of lower concen- 
tration (see §@). These properties are akin to those of 
LSB, gas-rich dlrr galaxies which reside in the outskirts 
of the LG, a fact that highlights tidal stirring as a plausi- 
ble causal mechanism for the origin of the morphology- 
density relation. 

7. The robustness of dwarf galaxies to tides and mass loss 
is increased significantly for those with more massive 
and more compact disks, as well as for dwarfs embed- 
ded in halos of higher concentration (see § 3J. This 
enhanced resilience to tidal effects has important con- 
sequences for the missing satellites problem as well as 
for determining the radial distribution of satellites in- 
side host halos (e.g., Diemand et al. 2 0041) . 
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8. The products of tidal stirring satisfy the relation V max = 
V^cr*, where cr* is the one-dimensional, central stel- 
lar velocity dispersion and V max is the maximum halo 
circular velocity (see § 15. It . Such a small conversion 
factor between er* and V max , formally valid for a tracer 
stellar population with an isotropic velocity dispersion 
tensor, is in agreement with those original ly adopted to 
formulate the missing sa tellites problem ( Moo re et al.l 
ll999HKlvpin et al.|[l999b . 

9. The mass-to-light ratios, M/L, of the orbiting disky 
dwarfs galaxies decrease monotonically with time 
as the extended DM halos are preferentially tidally 
stripped. In some cases, however, the M/L ratios start 
to increase later in the evolution when stellar mass loss 
becomes more effective (see § |5.3t . These cases are as- 
sociated with either enhanced tidal mass loss (R2 and 
R4) or more extended initial stellar distributions (R13) 
and demand that the dwarfs are stripped down to the 
scales where the alignment between the angular mo- 
menta of the stars and the orbital angular momenta of 
the dwarfs has an important effect on stellar stripping. 
This mechanism causes only a moderate increase of the 
M/L ratio and it may thus not be able to account for the 
extreme DM content in some of the classic dSphs (e.g., 
Draco or Ursa Minor). Producing such systems likely 
requires hydrodynamical processes. 

10. Distant dSphs in the LG, such as Leo I, Tucana, and Ce- 
tus, which are likely moving on very wide orbits, should 
have only been partially stirred by their hosts, assuming 



that their properties originate from tidal stirring. As a 
result, these remote dwarfs should exhibit higher values 
of Vrot/ c* compared to those of dSphs located closer to 
the primary galaxies. Future conclusive measurements 
of rotation in these systems will serve to validate (or 
falsify) this prediction. 
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